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ABSTRACT 

Context. Many aspects of massive star (^10 M Q ) formation are still unclear. In particular, the outflow properties at close distance 
(100-1000 AU) from a "Massive Young Stellar Object" (MYSO) are not yet well established. 

Aims. This work presents a detailed study of the gas kinematics towards the "Hot Molecular Core" (HMC) G31.41+0.31. 
Methods. To study the HMC 3-D kinematics at milli-arcsecond angular resolution, we have performed multi-epoch VLBI observations 
of the H 2 22 GHz and CH 3 OH 6.7 GHz masers, and single-epoch VLBI of the OH 1.6 GHz masers. 

Results. Water masers present a symmetric spatial distribution with respect to the HMC center, where two nearby (Ct.'2 apart), 
compact, VLA sources (labeled "A" and "B") are previously detected. The spatial distribution of a first group of water masers, named 
"Jl", is well fit with an elliptical profile, and the maser proper motions mainly diverge from the ellipse center, with average speed of 
36 km s -1 . These findings strongly suggest that the "Jl" water maser group traces the heads of a young (dynamical time of 1.3 10 3 yr), 
powerful (momentum rate of = 0.2 M Q yr~' km s~'), collimated (semi-opening angle ^ 10°) jet emerging from a MYSO located 
close (within «0'.'15) to the VLA source "B". Most of the water features not belonging to "Jl" present an elongated («2" in size), 
NE-SW oriented (PA«70 o ), S-shape distribution, which we denote with the label "J2". The elongated distribution of the "J2" group 
and the direction of motion, approximately parallel to the direction of elongation, of most "J2" water masers suggests the presence 
of another collimated outflow, emitted from a MYSO placed near the VLA source "A". The proper motions of the CH3OH 6.7 GHz 
masers, mostly diverging from the HMC center, also witness the expansion of the HMC gas driven by the "Jl" and "J2" jets. The 
orientation (PA«70°) of the "J2" jet agrees well with that (PA = 68°) of the well-defined Vlsr gradient across the HMC revealed 
by previous interferometric, thermal line observations. Furthermore, the "J2" jet is powerful enough to sustain the large momentum 
rate, 0.3 M yr -1 km s" 1 , estimated from the interferometric, molecular line data in the assumption that the Vlsr gradient represents 
a collimated outflow. These two facts lead us to favour the interpretation of the V L sr gradient across the G3 1.4 1+0.31 HMC in terms 
of a compact and collimated outflow. 

Key words. ISM: individual objects (G31.41+0.31) - ISM: kinematics and dynamics - masers - techniques: interferometric 



1. Introduction 

Understanding the formation of massive (OB-type) stars in- 
volves several theoretical and observational issues. From a theo- 
retical point of view, models of massive star-formation are com- 
plicated by the need to consider the effects of the intense stellar 
radiation of a MYSO, which, by heating, ionizing and exerting 
pressure on the circumstellar gas, strongly influences the pro- 
cess of mass accretion. For a spherically symmetric infall, the 
strong radiation pressure may halt and even reverse the accret- 



ing flow and thus stop further growth of the stellar mass ( Yorke 
2002|l. Recent three-dime nsional radiation-hydrodynamic mod- 



els (jKrumholz et al.|2009||"Kuiperet al.|2010||201 l[|Cunning ham 
et al.||2011| i indicate that even the most massive stars can be 
formed via accretion through a circumstellar disk, since the ra- 
diation pressure is decreased by conveying stellar photons along 
the disk axis and the ram pressure of the accreting gas increased 
by concentrating the accretion flow through the small disk solid 
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angle. However, magnetohydrodynamic calculations (Seifried 
et al.||20PT i suggest that the formation of a Keplerian disk can 
be suppressed by strong magnetic fields, since magnetic brak- 
ing can efficiently remove the excess angular momentum of the 
accreting gas, lowering the centrifugal support against gravity. 
Another class of models ("competitive accretion", |Bonnell &| 
|Bate| ( |2006| l), based on the observational evidence that massive 
stars are mainly found at the center of stellar clusters, predict that 
massive protostars accrete not only the gas of the surrounding 
core but also more diffuse gas from the whole molecular clump 
harbouring the cluster. According to this theory, accretion onto 
the massive protostar is driven by the gravitational potential well 
which funnels the gas down to the cluster center. In this scenario, 
circumstellar disks might be transient features prone to fragmen- 
tation and dispersion owing to strong interactions with cluster 
members. 

The most advanced simulations of massive star forma- 
tion model both the accretion and the ejection of mass self- 
consistently. The calculations by Cunningham et al. (2011 1 are 
the first to consider also the feedback effects of protostellar 
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outflows, as well as protostellar radiative heating and radiation 
pressure on the infalling, dusty gas. Previous works have con- 
sidered how ionization and strong radiation pressure can affect 
the mass ejection, suggesting that the resulting increase in the 
plasma pressure and turbulence could contribute to the decolli- 
mation of molecular outflows from massive OB-type protostars 
( |Konigl|[T999l |Yorke & Sonnhalter||2002l |Vaidya et al.|[20TT] l. 
Since massive stars very often belong to multiple stellar systems, 
strong member interactions could in principle influence the out- 
flow properties. In particular the tidal interaction between the 
disk and a companion star in a non-coplanar orbit has been pro- 



posed as a viable mechanism for jet precession (Terquem et al. 
[19991 |Bate et al.|2000] i. 

From an observational point of view, the two main difficul- 
ties to face for a study of massive star formation are 1) the rela- 
tively large distance (typically several kiloparsec) of the targets 
and 2) the large extinction at visible and (often) near-infrared 
wavelengths. The second point precludes direct observation of 
the (proto)star, whose properties can be inferred only through 
the study of the radio to mid-infrared radiation from the circum- 
stellar gas reprocessing the stellar photons. The first problem re- 
quires angular resolutions <0.' 1 to observe the scales (less than 
a few hundreds AU) relevant to the study of mass accretion and 
ejection. Prior to the advent of the Atacama Large Millimeter 
Array (ALMA), the angular resolution (> a few 07 1) of millime- 
ter interferometers (PdBI, SMA) was in most cases inadequate 
both to resolve massive accretion disks (with model-predicted 
sizes of a few hundreds AU) and to identify the protostellar out- 
flow emitted from a "single" MYSO system (hidden in the com- 
plex emission pattern resulting from the interaction of the cluster 
outflows). 

The limited angular resolution can explain the failure to de- 
tect disks in O-type stars, and the small number of bona-fide 
detections obtained for B-type stars ( |Cesaroni~e t al. 2007|. In 



association with the most luminous MYSOs one finds only huge 
($0.1 pc), massive (a few 100 M ) cores, with velocity gra- 
dients suggesting rotation. These objects, named "toroids", are 
likely non-equilibrium structures, because the ratio between the 
accretion time scale and the rotation period is much smaller 
in these toroids than in accretion disks (Cesar oni et al.| |2007; 

"" 201 The general properties of molecular out- 



Beltran et al. 



flows in massive star-forming regions are presently known 
mainly through single-dish surveys, with ^10" angular resolu- 



tion (Beuther et al. 


2002 Wu et al. 2005 Zhang et al. 2005 


Lopez-Sepulcre et a 


.|2009 1, of typical outflow tracers (such as 



CO and its isotopomer J CO). These studies have emphasized a 
good correlation between the momentum rate of the outflow and 
the MYSO bolometric luminosity. On the other hand, Very Large 
Array (VLA) sub-arcsecond observations have revealed a few 
thermal jets (HH 80-8 1 ( |Marti et al.|1998) ; IRAS 16547-424 7 
( [Rodriguez etal.|2008) ; IRAS 20126+4104 ( |Hofneretal.|2007) ) 
associated with intermediate- and high-mass (proto)stars. It is 
interesting that these (proto)stars, although being of moder- 
ate luminosity (10 4 — 10 s L Q ), eject very powerful, collimated 
thermal jets and molecular outflows, whose momentum rates 
(~0.1 M Q yr~' km s -1 ) are close to the upper limit measured 
in the single-dish surveys of molecular outflows. These results 
call for high-angular resolution ($0/1) observations of a larger 
sample of outflows from MYSOs to ascertain if the most com- 
pact and youngest flows (not observed in the single-dish surveys) 
are intrinsically more powerful and collimated. If this were con- 
firmed, the outflow evolution for high-mass (proto)stars would 
proceed similarly to the low-mass case, where outflows from 



Class objects are more intense and collimated than for Class I 
sources (Bontemps et al. 1 1 996 > . 

Since several years, we are conducting a campaign of multi- 
epoch Very Long Baseline Interferometry (VLBI) observations 
of the strongest interstellar molecular masers (OH at 1 .6 GHz, 
CH 3 OH at 6.7 GHz, and H 2 at 2 2 GHz) towards a small sample 
of massive star-forming regions (Moscadelli et al. 2 007| 201 1 1 



|Sanna et al.||2010a[ |Goddi et al.f|2011||Li et al.||2012) . Multi- 
epoch VLBI observations, achieving angular resolutions of 1- 
10 mas, permit to measure the proper motions of the maser emis- 
sion centers, and, knowing the source distance and the maser 
Vlsr (from Doppler effect), one can reconstruct the 3-D kine- 
matics of the masing gas. Our results demonstrate the ability 
of maser emission to trace kinematic structures close to (within 
tens/hundreds AU from) the MYSO, revealing the presence of 
fast wide-angle and/or collimated outflows, traced by the H2O 
masers, and of rotation and infall, indicated by the CH3OH 
masers. This work reports on our recent VLBI observations of 
the molecular masers associated with the HMC G31.41+0.31, 
whose gas kinematics has been investigated in the last years us- 
ing the best available cm and mm interferometers (VLA, SMA, 
PdBI). 

G3 1.41 +0.31 is located at a kinematic distance o f 7.9 kpc, 
and its high bolometric luminosity of 3X10 3 L Q (Cesaroni 



et al. 1994b 1 hints at the presence of embedded high-mass 
(proto)stars. The HMC in G31.41+0.31 (at a systemic LSR ve- 
locity, V sys = 97.0 km s~') was originally imaged with the VLA 
in the high-excitation (4,4) inversion transition of ammonia by 



Cesaroni et al. (1994ai, and successively mapped by Beltran 



etaL[(2004) and|CesaronietaI^([20Tl), using the PdBI and SMA 



interferometers, respectively, in the (12-11) rotational transitions 
of methyl cyanide (CH3CN). These latter observations, achiev- 
ing an angular resolution of 0/8, reveal a well-defined veloc- 
ity gradient across the HMC, oriented in the NE-SW direction. 
Considering the large size (»0.1 pc) and mass (several 100 M Q ) 
of this structure, Beltran et al. (2004)l and |Cesaroni et al.] ( |201 \\ 
interpret the velocity gradient in terms of a toroid in pseudo- 
Keplerian rotation around a stellar cluster of ^10 3 M Q , tightly 
packed at the center of the HMC. In contrast to this interpreta- 
tion, [Arayaetal]([2008 ) explain the velocity gradient in terms of 
a compact outflow on the basis of the NE-SW bipolar morphol- 
ogy of thermal methanol (at 44 GHz), observed with the VLA at 
an angular resolution of a;0/5. In an attempt to decipher the gas 
kinematics inside the G31.41+0.31 HMC, we have observed the 
OH 1.6 GHz, CH3OH 6.7 GHz, and H 2 22 GHz masers asso- 
ciated with the core, using the VLBI. The observational details 
are given in Sect. [2] while the results are illustrated in Sect. [3] 
and discussed in Sects. [4] and [5] The conclusions are drawn in 
Sect. [6] 

2. VLBI Maser Observations and Data Reduction 

We conducted VLBI observations of the H 2 22 GHz 
and CH3OH 6.7 GHz masers (at four epochs) and of the 
OH 1.665/1.667 GHz masers (at a single epoch) toward 
G3 1.41 +0.31. Table [T] lists observing epochs and beams, and 
reports the derived maser positions. In order to determine the 
maser absolute positions, we performed phase-referencing ob- 
servations in fast switching mode between the maser source 
and the calibrator J1834-0301. This calibrator has an angu- 
lar offset from the maser source of 3?8 and belongs to the 
list of sources of the International Celestial Reference Frame 2 
(ICRF2). Its absolute position is known with a precision better 
than +1.0 mas, and its flux measured at S and X bands is 250 and 
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280 mJy beam , respectively. Four fringe finders, J1642+3948, 
J1751+0939, J1800+3848 and J2253+1608, were observed for 
bandpass, single-band delay, and instrumental phase-offset cali- 
bration. 

Data were reduced with AIPS following the VLBI spectral 
line procedures. For a description of the general data calibration 
and the criteria used to identify individual masing clouds, derive 
their parameters (position, intensity, flux and size), and measure 
their (relative and absolute) proper motions, we refer to the re- 
cent paper on VLBI observations of H2O and CH3OH masers 



by |Sanna et al.| ( |2010a| ). For each VLBI epoch, except for the 

, H z O and CH 3 OH 



fourth epoch on CH3OH masers (see Sect |2.3 
maser absolute positions have been derived by using two ways, 
direct phase-referencing and inverse phase-referencing, and the 
two procedures always gave consistent results. For the weak 
OH masers, we only use direct phase-referencing to determine 
their absolute positions, because inverse phase-referencing is 
hindered by the insufficient number of maser (phase calibra- 
tion) solutions. The derived absolute proper motions have been 
corrected for the apparent proper motion due to the Earth's or- 
bit around the Sun (parallax), the Solar Motion and the dif- 
ferential Galactic Rotation between our LSR and that of the 
maser source. We have adopted a flat Galax y rotation curve 
(Rn = 8.3 ± 0.23 kpc, @ = 239 + 7 km s _1 ) (jBrunthaler et al 



2011 



, and the Solar motion ( U = ll.ljffiS 



) byjSchonrich et al. (|2010|l, who have re- 



V = 12.24^-^, and 



W = 7.25t°;K km s" 

cently revised the Hipparcos satellite results. We antecipate that 
the absolute proper motions of water and methanol masers are 
dominated by a large systematic effect due to the uncertainties 
in the source distance and Solar and Galactic motion, and/or by 
a peculiar velocity of the maser source with respect to its LSR 
reference system. Being severely affected by this systematic de- 
viation, the absolute proper motions cannot be used to describe 
the maser kinematics around the MYSO, to which purpose we 
employ the relative proper motions. 



2.1. VLBA: 22.2 GHz H 2 Masers 

We observed G31.41+0.31 (tracking center: R.A.(J2000) = 
18 h 47 m 34.35 s and Dec.(J2000) = -01°12'45.6") with the 
NRACQ Very Long Baseline Array (VLBA) in the 6 i6 - 5 23 
H2O transition (rest frequency 22.235079 GHz). The obser- 
vations (program code: BM286) consisted of 4 epochs: 2008 
November 23, 2009 January 18, 2009 March 18, and 2009 May 

( fT999l >, 



18. Following Reipurth et al. ( |1999[ ), to reduce systematic er- 
rors owing to correlator atmospheric mismodeling, we observed 
three "geodetic" blocks, each lasting about 40 min, before the 
start, in the middle, and after the end of our phase-reference 
observations. The "geodetic -type" data were taken in left circu- 
lar polarization with eight 4 MHz bands spanning a bandwidth 
of 480 MHz. The fast-switching phase-referenced observations 
consisted of two «3 hour blocks, during which we recorded the 
dual circular polarization through a 16 MHz bandwidth centered 
at Vlsr = 96 km s _1 . The data were processed with the VLBA 
FX correlator in Socorro (New Mexico) using an averaging time 
of 1 s and 1024 spectral channels. 

Images were produced with natural weighting, cleaned and 
restored with the naturally-weighted beam, having a Full Width 
at Half Maximum (FWHM) size of about 1.6 mas x 1.0 mas 
at a position angle (PA) of 6° (East of North), showing small 



variations from epoch to epoch (see Table[T|. The interferometer 
instantaneous field of view was limited to ~2'.'5. In each observ- 
ing epoch the on-source integration time was ~2.5 h, resulting 
in an effective rms noise level on the channel maps (lcr) in the 
range 6-10 mJy beam -1 . The spectral resolution was 0.21 km 



2.2. VLBA: 1 .665/7 GHz OH Masers 

We observed G31.41+0.31 (tracking center: R.A.(J2000) = 
18 h 47 m 34.3 s and Dec.(J2000) = -01°12'46.0") with the VLBA 
in the "main" lines (rest frequencies of 1.665401 GHz and 
1.667359 GHz) of the hyperfine splitting of the 2 Y1 3/2 J = 3/2 
OH level on September 26 2008 (program code: BM286E). 
During a run of ~8 h, we alternated the observations of either 
"main" lines, recording the dual circular polarization with one 
bandwidth of 1 MHz (centered at the maser Vlsr of 97.0 km 
s _1 ) and three 4 MHz bandwidths. The 4 MHz bandwidths were 
used to increase the SNR of the weak L-band signal of the con- 
tinuum calibrator. The data were processed with the VLBA FX 
correlator in two correlation passes using 1024 and 128 spectral 
channels for the 1 MHz and 4 MHz bands, respectively. In each 
correlator pass, the data averaging time was 2 s. 

Images were produced with natural weighting, cleaned and 
restored with the naturally-weighted beam, having a FWHM size 
of 19 mas x 1 1 mas at a PA of -3°, and 24 mas X 15 mas at a PA 
of 35°, for the 1665 and 1667 MHz masers, respectively. The in- 
terferometer instantaneous field of view was limited to 7'.' 2. The 
on-source integration time was ~1.2 h, resulting in an effective 
rms noise level in each velocity channel of ~17 mJy beam -1 . 
The spectral resolution attained across the maser 1 MHz band 
was 0.18 km s _1 . 



2.3. EVN: 6.7 GHz CH 3 OH Masers 

We observed G31.41+0.31 (tracking center: R.A.(J2000) = 
18 h 47 m 34.33 s and Dec.fJ2000) = -01°12'47.0") with the 
European VLBI Network (EVN) in the 5i - 6 A + CH 3 OH 
transition (rest frequency 6.668519 GHz). Data were taken at 
4 epochs (program code: EM071) separated by ~1 year: 2009 
March 15, 2010 March 11, 2011 March 2, and 2012 February 
24. In the first two epochs, the antennae involved in the observa- 
tions were Cambridge, Jodrell2, Effelsberg, Onsala, Medicina, 
Noto, Tomn, Westerbork and Yebes. Because of technical prob- 
lems, the Noto antenna could not take part to the observations 
of the third and fourth epochs. The Knockin antenna replaced 
the Cambridge antenna in the third epoch, and the Cambridge 
antenna did not observe the fourth epoch. During a run of ~8 h, 
we recorded the dual circular polarization through eight adjacent 
bandwidths of 2 MHz, one of them centered at the maser Vlsr of 
99.0 km s _1 . The eight 2 MHz bandwidths were used to increase 
the SNR of the weak continuum calibrator. The data were pro- 
cessed with the MKIV correlator at the Joint Institute for VLBI 
in Europe (JIVE-Dwingeloo, the Netherlands) in two correla- 
tion passes, using 1024 and 128 spectral channels to correlate 
the maser 2 MHz bandwidth and the whole set of eight 2 MHz 
bandwidths, respectively. In each correlator pass, the data aver- 
aging time was 1 s. 

In addiction to the absence of Cambridge and Noto, the 
fourth observing epoch was also severely affected by other prob- 



1 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with the 
National Science Foundation. 



The European VLBI Network is a joint facility of European, 
Chinese and other radio astronomy institutes founded by their national 
research councils. 
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lems: the Jodrell2 receiver was warm (T sys » 250 K); Medicina 
suffered a computer problem, after which the fringe ampli- 
tudes were much reduced; because of strong winds, Onsala was 
stowed half of the observing time. The good data were not suffi- 
cient to reliably measure the maser absolute position, which has 
been merely estimated from a linear fit of the positions of the 
first three epochs (see Table[T]l. 

Reflecting the changes in the antennae taking part in the ob- 
servations, the FWHM size of the naturally-weighted beam var- 
ied over the four epochs from about 8 mas x 5 mas (PA = -17°) 
to about 12 mas x 3.5 mas (PA = -37°) (see Table[T}. At the first 
two epochs, the maser images have been cleaned and restored us- 
ing the naturally-weighted beam. At the third and fourth epoch, 
to prevent apparent changes in the maser structure owing to the 
variation of the naturally-weighted beam, we have restored the 
maser images using a beam, 8.3 mas x 4.8 mas at PA = -17°, 
average of the naturally-weighted beams of the first and second 
epochs. The interferometer instantaneous field of view was lim- 
ited to ~33". Using an on-source integration time of ~2.3 h, the 
effective rms noise level on the channel maps (lcr) varied in the 
range 5-13 mJy beam -1 . The spectral resolution attained across 
the maser 2 MHz band was 0.09 km s . 



3. Results 

3.1. 22 GHz Water Masers 

The lower panel of Fig. [T] shows absolute positions and LSR 
velocities of the detected water maser features, whose param- 
eters (derived by fitting a 2-D elliptical Gaussian to the maser 
intensity distribution) are listed in Table |4] Over four epochs, 
we detected 173 distinct maser features, whose intensity ranges 
from the detection threshold of »50 mJy beam 1 up to 
~100 Jy beam -1 , and with Vlsr varying from 84 km s _I to 
128 km s . Interestingly, the distribution of the water masers 
is approximately symmetric with respect to the position of the 
two nearby, compact VLA continuum sources "A" and "B", ob- 
served by Cesaroni et al. (|2010|> roughly at the center of the 



G31.41+0.31 HMC. At comparable offsets to the North and 
South of the continuum sources, one observes two similar clus- 
ters of water masers (enclosed in the white dotted rectangles in 
Fig. [T}, both clusters having the shape of acute bows, with the 
bow tip pointing to the North (South) for the northern (southern) 
cluster. In the following, we will refer to these two maser clusters 
with the name "Jl", and will denote the northern and southern 
"Jl" clusters with "Jl-N" and "Jl-S", respectively. Most of the 
water features not belonging to "Jl" present a NE-SW oriented, 
S-shape distribution (inside the white dashed polygon of Fig.[TJ, 
which consists of a NE-SW extended («2"in size, PA»70°) strip 
of features, and two smaller (w0'/5 in size), SE-NW oriented 
lines of masers at the ends of the NE-SW strip. In the follow- 
ing, we refer to this S-shape pattern of maser features as "J2", 
and use the names "J2-SW" and "J2-NE" to indicate the water 
maser clusters at the southwestern and northeastern ends of the 
"J2" distribution, respectively. We note that while the Vlsr of 
the "J2" masers are relatively close (within ±10 km s -1 ) to V sys , 
"Jl" includes the most blue- and red-shifted water features (up 
to 30 km s _1 away from V sys ). 

We calculated the geometric center (hereafter "center of mo- 
tion", identified with label #0 in Table |4) of the features with a 
stable spatial and spectral structure, persisting over the four ob- 
serving epochs, and refer our measurements of proper motions 
to this point. The upper panel of Fig. [2] presents the proper mo- 
tions of the water masers relative to their "center of motion". The 



magnitude of the relative proper motions (measured for 70 fea- 
tures and listed in Tabled ranges from »10kms 'to » 120 km 
s _1 , with a mean value of 34.0 km s~', and a mean error of 
7.0 km s -1 . The proper motions of the water masers belonging 
to the "Jl-N" and "Jl-S" clusters indicate expansion, since most 
of the proper motion vectors diverge from the HMC center and 
are oriented approximately perpendicular to the pattern outlined 
by the water masers. Concerning the velocity distribution of the 
"J2" water masers, the few measured proper motions for the fea- 
tures found in the NE-SW strip suggest that the masers are re- 
ceding from the HMC center moving along directions at close 
angle with the strip orientation. The water features in the "J2- 
S W" cluster appear to move consistently towards SW, parallel to 
the strip orientation. In contrast, the proper motions of the "J2- 
NE" maser features present a more scattered distribution, with 
the direction of motion rotating from NE to NW with the maser 
position varying from SE to NW. Finally, we note that the water 
maser features not belonging to "Jl" and "J2" groups, clustered 
to the SE of the whole maser distribution, present motions to the 
North. 

To check the reliability of the relative, water maser proper 
motions we follow the same procedure applied in previous stud- 
ies. Assuming that in G3 1.41 +0.31 (as verified in several sources 
dMoscadelli et~aT1[20T0l |20TT] |Sanna et"aTT|2010a|b) ) methanol 
masers move (with respect to the star) significantly slower than 
the water masers, we can use the absolute proper motion of 
the methanol masers to correct the absolute proper motion of 
the water masers. The result should be a good approximation 
of the water maser velocities with respect to the reference sys- 
tem comoving with the star. The lower panel of Fig. [2] shows 
the absolute water maser proper motions after subtracting the 
average absolute, methanol maser proper motion (evaluated in 
Sect.. [373) . Note that, within the errors, the methanol-corrected 
absolute proper motions of the water masers are consistent with 
the relative proper motions calculated with respect to the "cen- 
ter of motion" (upper panel of Fig. [2|. This result supports our 
assumption that the "center of motion" can be a suitable refer- 
ence system to calculate water maser velocities. In the following 
discussion, we use the relative velocities of the water masers to 
describe the gas kinematics close to the MYSO. 

The methanol-corrected absolute proper motions reported in 
the lower panel of Fig. [2] should not be confused with the ab- 
solute proper motions of the water maser features, derived from 
the measured absolute proper motion (corrected for the apparent 
motion) of feature #1 (see Table[T|. The absolute proper motions 
differ from the relative proper motions presented in the upper 
panel of Fig. |2]by a vector of amplitude ^60 km s~' pointing 
to NE, which seems to be the dominant component of all proper 
motions and reflect the large systematic deviation affecting the 
measurement of the absolute proper motions. 

3.2. 1 .6 GHz Hydroxy I Masers 

The lower and upper panels of Fig. [3] show absolute positions 
and LSR velocities of the detected 1665 MHz and 1667 MHz 
OH maser features, respectively, whose parameters are listed 
in Tables [2] and [3] The number of detected features is sig- 
nificantly higher at 1667 MHz (52) than at 1665 MHz (17), 
consistent with the fact that the average OH maser intensity 
is higher at 1667 MHz (0.9 Jy beam" 1 ) than at 1665 MHz 
(0.7 Jy beam -1 ). For both OH maser transitions, the detec- 
tions are similarly distributed between RCP and LCP compo- 
nents, which also show comparable intensities. The overall OH 
maser Vlsr interval extends from 89 km s" 1 to 111 km s"" 1 . The 



4 



L. Moscadelli et al.: A Double- Jet System in the G3 1.4 1+0.31 Hot Molecular Core. 

Vlsr (km s" 1 ) 




J i i i i I i i i i I i i i i L 

18 h 47 m 34740 34735 34?30 34 s 25 

RIGHT ASCENSION (J2000) 

Fig. 1. Lower and upper panels show the absolute position and LSR velocity of the 22 GHz water (filled triangles) and 6.7 GHz methanol (dots) 
maser features, respectively. Colors denote the maser Vlsr according to the color-velocity conversion code shown on the right-hand side of the 
upper panel. Symbol area is proportional to the logarithm of the feature intensity. In both panels, the grey-scale map reproduces the integrated 
emission of the CH 3 CN (12-11) K=4 line (at »221 GHz) observed with the SMA (beam k078) by |Cesaroni et al.| < [20TT| >. Contour levels are from 
0.15 Jy beam 1 to 2.3 Jy beam - ' by steps of 0.27 Jy beam -1 . The white big dots denote the position of the two VLA compact sources, named 
"A" and "B", detected by Cesa roni et al.| ( (2010j >. In both panels, the two dotted rectangles and the S-shape dashed polygon encompass the area 
where the "Jl" and "J2" water maser groups distribute, respectively (see Sect |3.1[ l. The northern and southern "Jl" clusters are named "Jl-N" and 
"Jl-S", respectively. The water maser clusters at the southwestern and northeastern ends of the "J2" distribution are named "J2-SW" and "J2-NE", 
respectively. In the upper panel, the dotted ellipse surrounds a cluster of 6.7 GHz methanol masers, identified with the name "F", more detached 
from the area of water maser detection. c 
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Fig. 2. Upper and lower panels show the relative (referred to the "center of motion") and (methanol-corrected) absolute proper motions of the 
22 GHz water masers, respectively. Colored vectors indicate the direction and the amplitude of the proper motions. The black vector in the bottom 
left corner of the panels indicates the amplitude scale of proper motions in kilometer per second. Dotted vectors are used for the most uncertain 
proper motions. Colors, symbols, the grey-scale map, the white big dots, the dotted rectangles, and the S-shape dashed polygon, have the same 
meaning as in Fig.^ 
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Table 1. Maser positions, brightnesses and corresponding synthesized beams. 



Maser 


Feature 


Epoch of 


R.A.(J2000) 


Dec.(J2000) 


Ipeak 


Major x Minor, PA 




Number 


Observation 




(° ' ") 


(Jy beam -1 ) 


(mas x mas, deg) 


H z O 


1 


2008 Nov 23 


18 47 34.31259 


-01 12 46.6384 


115 


1.6x0.8,-7 






2009 Jan 18 


18 47 34.31257 


-01 12 46.6394 


144 


1.6 x 1.1 ,-1 






2009 Mar 18 


18 47 34.31256 


-01 12 46.6402 


75 


1.4 x 1.0 , 13 






2009 May 18 


18 47 34.31254 


-01 12 46.6412 


40 


1.8 x 1.1 ,20 


CH 3 OH 


1 


2009 Mar 15 


18 47 34.28652 


-01 12 45.5777 


8.3 


7.8x5.0,-17 






2010 Mar 11 


18 47 34.28636 


-01 12 45.5837 


9.0 


8.8x4.6,-18 






2011 Mar 2 


18 47 34.28625 


-01 12 45.5879 


10.9° 


11.8 x 3.3 , -36 






2012 Feb 24 


18 47 34.28610 


-01 12 45.5934 


10.1" 


11.6x3.3,-38 


OH- 1665 


1 LCP 


2008 Sep 26 


18 47 34.27583 


-01 12 46.3282 


4.0 


19 x 11 , -3 


OH- 1667 


1 LCP 


2008 Sep 26 


18 47 34.27574 


-01 12 46.3287 


5.0 


24 x 15 , 35 



Notes. Column 1 indicates the maser species; column 2 gives the label number (and, for the OH masers, polarization) of the phase-reference feature, 
as reported in Tables [2] [3] [4] and [5] column 3 lists the observing date; columns 4 and 5 report the absolute position in Equatorial coordinates; 
column 6 gives the peak intensity; column 7 reports the (major and minor) FWHM size and PA of the naturally-weighted beam. The PA of the 
beam is defined East of North. The accuracy in the derived maser absolute positions is ~ 1 mas in both Equatorial coordinates. 
<fl) Peak intensity from the image restored using a beam with FWHM size of 8.3 mas x 4.8 mas at PA = -17°. 



1665 MHz and 1667 MHz OH maser spatial distributions are 
similar, both extending more along the East- West (3"-3'.'5) than 
along the North-South (~ 1") direction. In particular, we note that 
the spatial distribution of both the 1665 MHz and 1667 MHz OH 
masers is not symmetrical with respect to the two VLA contin- 
uum sources found at the center of the HMC. Both maser transi- 
tions populate preferentially the region to the West of the VLA 
sources, spreading more or less uniformly up to a maximum off- 
set of *2" . Instead, to the East of the VLA sources, the observed 
OH maser features concentrate at relatively small (<075) offsets 
and very few features are found at large (~T.'5) separation from 
the HMC center. The OH masers present also a regular change in 
Vlsr with the position, with Vlsr increasing from West to East. 
We discuss further the properties of the OH maser spatial and 



velocity field in Sect. 4.3 



3.3. 6.7 GHz Methanol Masers 

The upper panel of Fig. [T] shows absolute positions and LSR 
velocities of the detected methanol maser features, whose pa- 
rameters are listed in Table [5] Over four epochs, we detected 
85 distinct maser features, whose intensity ranges from the de- 
tection threshold of »40 mJy beam -1 up to »T0 Jy beam" 1 , 
and with Vlsr varying from 90 km s -1 to 107 km s _1 . 6.7 GHz 
masers distribute at distances SSI" from the two compact VLA 
continuum sources (at the HMC center), mainly clustering to the 
SE, SW and NW of the continuum positions, with no 6.7 GHz 
features detected to the East of the VLA sources. Comparing the 
methanol and water maser spatial distributions, we find two clus- 
ters of 6.7 GHz features (red- and blue-shifted and placed inside 
the northern and the southern dotted rectangles of Fig. [T] respec- 
tively) which have positions and Vlsr similar to 22 GHz features 
belonging to the North-South oriented, "Jl" group. A resem- 
blance in positions and Vlsr between the two maser species is 
also observed to the SW of the VLA continuum sources, where 
methanol masers are spread nearby the "J2-SW" water maser 
group. We find, however, also a small cluster of 6.7 GHz features 
to the NW of the continuum sources (inside a dotted ellipse in 
Fig. [T}, detached from the area of water maser detection. In the 
following, we will refer to this 6.7 GHz maser cluster with the 
name "F". 

We calculated the geometric center (hereafter "center of mo- 
tion", identified with label #0 in Table [5]) of the features with 



a stable spatial and spectral structure, persisting over the 4 ob- 
serving epochs, and refer our measurements of proper motions to 
this point. The upper panel of Fig. [4]presents the proper motions 
of the methanol masers relative to their "center of motion". The 
magnitude of the relative_proper motions (measured for 40 fea- 
tures and listed in Table [5} ranges from «5 km s -1 to «30 km 
s , with a mean value of 13.0 km s _1 , and a mean error of 
4.0 km s 1 . It is interesting to compare the proper motions of the 
water and methanol masers. The two clusters of 6.7 GHz features 
in proximity of the "Jl" H2O maser structure (inside the dotted 
rectangles of Fig. [4] upper panel) have proper motions similarly 
oriented as the nearby water masers, with CH3OH features in 
correspondence of the northern (southern) "Jl" lobe, moving to 
NW (SE) (see Fig [5] right panel). The "Jl" water masers, how- 
ever, move on average faster than the adjacent methanol masers. 
Proper motions of only «5 km s -1 (pointing to N-NE) are mea- 
sured for the group of 6.7 GHz masers offset 07 1-075 to the 
West of the the "Jl-N" water maser cluster. The methanol masers 
with positions intermediate between the "Jl-S" and the "J2-SW" 
water maser clusters move preferentially to South, while the mo- 
tion of the closeby, "J2" water masers is mainly oriented to SW. 
Finally, we note that the features in the 6.7 GHz maser cluster 
"F" (inside a dotted ellipse in Fig.|4] upper panel) present proper 
motions in opposite directions. 



The lower panel of Fig. [4] shows the absolute proper mo- 
tions of the 6.7 GHz CH3OH masers. These differ from the rel- 
ative proper motions (presented in the upper panel of Fig. |4]i 
by a vector, corresponding to the average absolute proper mo- 
tion of the 6.7 GHz masers, with components of 29 ± 4 km s -1 
along the right-ascension axis and 29 ± 7 km s -1 along the dec- 
lination. As it is clear from Fig. |4] the dispersion (»T0 km s -1 
in each Equatorial coordinate) of the 6.7 GHz maser, absolute 
proper motions from the mean is small, suggesting, as noted in 
Sect. - 



3.1 



that a large systematic effect dominates our measure- 
ment of absolute proper motions. Since for the 6.7 GHz CH3OH 
masers, this systematic deviation (with amplitude of »41 km 
s -1 ) overwhelms the average amplitude (13 km s -1 ) of the rela- 
tive proper motions, its value can be adequately estimated as the 
average of the 6.7 GHz maser, absolute proper motions. As al- 
ready explained in Sect. 



3.1 



the absolute proper motions of the 
22 GHz H2O masers presented in the lower panel of Fig|2]have 
been corrected for this systematic effect. 



7 



L. Moscadelli et al.: A Double- Jet System in the G3 1.4 1+0.31 Hot Molecular Core. 



Vlsr (km s" 1 ) 



O 
O 
O 
W 
•-» 

sz; 

o 
i— i 

E- 
< 

a 

o 
w 

Q 



CD 



w _ 



CD 



1^ 



A 1667 RCP 
□ 1667 LCP 




1 r~ 

2000 AU 



H 1 h 



~k 1665 RCP 
O 1665 LCP 



2000 AU 




I8 h 47 m 34 s 4 



34 s 3 



34 s 2 



RIGHT ASCENSION (J2000) 



84 

85.3 

86.6 

87.9 

89.2 

90.5 

91.8 

93.1 

94.4 

95.7 

100.4 

103.9 

107.3 

110.8 

114.2 

117.7 

121.1 

124.6 

128.3 



Fig. 3. Lower and upper panels show the absolute position and LSR velocity of the 1665 MHz and 1667 MHz OH masers, respectively. For the 
1665 MHz masers, the LCP and RCP features are indicated with empty pentagons and stars, respectively. Empty squares and empty triangles 
indicate the LCP and RCP features of the 1667 MHz masers, respectively. Symbol area is proportional to the logarithm of the feature intensity. 
The black dashed line in the lower panel shows the axis of projection of the OH maser positions to produce the position-velocity plot of Fig. [6] 
Colors, the grey-scale map, and the white big dots, have the same meaning as in Fig. [7] 



4. A Double-Jet System inside the G31.41+0.31 
HMC 

4. 1 . The "J1 " Water Maser Group 



In Sect. 3. 1 we noted the symmetrical spatial and velocity distri- 
butions of the water masers belonging to the "Jl-N" and "Jl-S" 
clusters. Figure [5] (left panel) shows that these distributions can 
be well fit with an ellipse. The fit parameters are: 1) the position 
on the sky of the ellipse center; 2) the major and minor axes; 
3) the PA (East from North) of the major axis. The best-fit posi- 
tion of the ellipse center is determined to be 07 14+(y.'03 South 
(PA = 177°+5°) of the VLA source "B" (i.e. the NE source; 



see Fig. [TJ, which is closer to the "Jl" center than the VLA 
source "A". The best-fit major and minor axes are l'.'44±0'.'l 
and 0724±0702, respectively, with PA = 8°±1°. 

An elliptical spatial distribution is expected if the masers 
traced a circular ring inclined with respect to the plane of the 
sky. However, assuming that the masers also move on the ring 
plane (as, for instance, in an expanding ring), this scenario is 
clearly to be ruled out on the basis of the comparison of the fit 
ellipse eccentricity and the 3-D maser velocities. The ratio be- 
tween the minor and the major axis of the best-fit ellipse to the 
"Jl" maser distribution gives an inclination angle of the plane of 
the putative maser ring with respect to the plane of the sky of 
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Fig. 4. Upper and lower panels show the relative (referred to the "center of motion") and absolute (corrected for the apparent motion, see Sect.|2^ 
proper motions of the 6.7 GHz methanol masers, respectively. Colored vectors indicate the direction and the amplitude of the proper motions. The 
black vector in the bottom right corner of the panels indicates the amplitude scale of proper motions in kilometer per second. Dotted vectors are 
used for the most uncertain proper motions. Colors, symbols, the grey-scale map, the white big dots, the dotted rectangles, the S-shape dashed 
polygon, and the dotted ellipse, have the same meaning as in Fig.[TJ 
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Fig. 5. Each of the three panels shows the same area of the sky around the "Jl" water maser group. Maser positions are always refered to the 
VLA source "B" (see text). In each panel, the dashed line shows the best-fit ellipse to the distribution of the "Jl" water masers. Left panel: Black 
triangles denote the positions of the "Jl" water masers. The dotted lines trace the major (labeled "P-axis") and minor ("T-axis") axis of the best-fit 
ellipse. Black dots give the positions of the VLA sources "A" and "B". Central and Right panels: Colored triangles and dots show the positions 
of the water {central panel) and methanol {right panel) masers, respectively. Symbol area is proportional to the logarithm of the feature intensity. 
Colors, the grey -scale map, the white big dots have the same meaning as in Fig. [I] Colored vectors indicate the direction and the amplitude of the 
maser proper motions. The black vector in the bottom right corner of the panels indicates the amplitude scale of proper motions in kilometer per 
second. Dotted vectors are used for the most uncertain proper motions. 



/ c i r = arccos(0724/r.'44) 80°, indicating that the maser circle 
should be seen close to edge-on. On the other hand, if masers 
move on the ring plane, the ring inclination angle can be also 
evaluated by the ratio of maser velocity components projected 
along the line-of-sight, V z , and along the ellipse minor axis, Vj. 
Averaging on all the (39) "Jl" features with measured proper 
motions, since in general V r » V z , the ring inclination angle 
as determined from the maser velocities is 2°±2°, implying that 
the maser motion should be seen face-on. That contrasts with the 
elongated maser spatial distribution, and we have to rule out the 
ring interpretation (with masers moving on the ring plane). 



If the H2O masers in "Jl" are accelerated by a bipolar jet di- 
rected approximately North-South, that could explain the elon- 
gated elliptical distribution of the masers and the pattern of 
proper motions mainly diverging from the ellipse center. The pu- 
tative location of the MYSO driving this jet likely coincides with 
the ellipse center. As previously noted, the best-fit position of the 
ellipse center lies only »0714 to the South of source "B" (see 
Fig. [5] left panel). Moreover, based on the cm-wavelength SED 



analysis, |Cesaroni et al. 



( |2010| l reason that the source "B" could 



be a thermal jet rather than a compact H II region. We conclude 
that source "B" could be excited by the same (proto)stellar jet 
driving the "Jl" water maser motions. 



One can estimate the momentum rate of the jet from the av- 
erage distance of the H2O masers from the MYSO and the av- 
erage maser velocity. Assuming a H2 pre-shock density «h 2 = 
10 8 cm" 3 , as predicted by excitation models of 22 GHz masers 
( Elitzur et al.|19 89|, the momentum rate in the water maser jet 
is given by the expression: 



P = 1.5X 10" 



■V^ R 2 m (n/47i) Moyr-'kms" 1 



(1) 



Here V10 is the average maser velocity in units of 10 km s" 1 , 
R\oo the average distance of water masers from the YSO in units 
of 100 AU, and O the solid angle of the jet. This expression has 
been obtained by multiplying the momentum rate per unit sur- 
face transferred to the ambient gas («h 2 m H 2 V 2 )> by D.R 2 , un- 
der the assumption that the jet is emitted from a source at a 
distance R from the masers within a beaming angle Q.. To es- 
timate the average distance and velocity of the water masers 
tracing the jet, we use all the (39) maser features with mea- 
sured proper motions belonging to the "Jl" group. The result- 
ing average distance is 0V47 (or »3700 AU) and the average 
speed is 36 km s~'. Using these values, the jet momentum rate 
is P a 27 (£1/ An) M Q yr~ 1 kms 4 . 

The ratio between the minor and the major axis of the best-fit 
ellipse can be used to estimate the jet semi-opening angle 0j a 
10°, which, assuming that the jet is symmetric about its axis, cor- 
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responds to a jet solid angle Q ^ 0.1 sterad. Using such a value 
for Q, the jet momentum rate results P =t 0.2 M Q yr~' km s _1 . 

Figures [T] and |2] show that, in both the "Jl-N" and "Jl-S" 
clusters, most of the water maser features are found projected 
close to the peaks of the bulk CH3CN (12-11) line emission. A 
possible interpretation is that the H2O masers are excited at the 
interface between the jet and the surrounding dense molecular 
core. The profile outlined by the water masers in the "Jl-N" and 
"Jl-S" clusters reminds that of the bow-shocks associated with 
Herbig-Haro objects (Hartig an et al.|20TT) , produced by proto- 
stellar jets when they impinge against the dense circumstellar 
material. This scenario suggests that the "Jl" jet is still confined 
inside a dense cocoon of dust and gas within a distance <1" 
from the MYSO, which is further indicated by the fact that, de- 
spite our sensitive observations, we have failed to reveal a molec- 
ular outflow along the "Jl" (North-South) direction up to angular 
scales of several 10"( |Cesaroni et aL|2011) . Using the average 
value of »20 km s for the jet-parallel component of maser ve- 
locities, and the jet major semi-axis of 077, one derives a short 
dynamical time scale of 1.3 10 3 yr. 

4.2. The "J2" Water Maser Group 

Looking at the "J2" water maser distribution (see Figs.[T]and[2|, 
the NE-SW elongated strip of water masers and the motion of 
the water features in the "J2-SW" cluster, parallel to the strip 
orientation, suggests the presence of another collimated outflow, 
at PA«70°. We speculate that the water masers in the strip mark 
an ionized jet, amplifying the free-free continuum on the back- 
ground: in this view the emission from source "A", which lies 
along the maser strip (see Fig. [TJ, could be tracing the jet close 
to the MYSO. However, the proposed scenario cannot readily ac- 
count for the direction of motion of water masers in the "J2-NE" 
group, where several maser features move to the North rather 
than parallel to the direction (PA»70 o ) of the putative jet. Note 
also that to the South (offset by of the "J2-NE" cluster, an- 
other cluster of water maser features (inside the dashed circle in 
Fig. [7] lower panel) with proper motions directed close to North 
is found. The motion of these features do not seem to be consis- 
tent with the jet interpretation for the "J2" masers. Section 5.2 
propones a qualitative explanation of the more scattered velocity 
distribution for the masers to the East of the VLA source "A", 
considering the dynamical interaction of the HMC gas with its 
environment. In the following we discuss the scenario that the 
"J2" masers trace a jet. 

Repeating the calculation performed for the "Jl" maser 
group in Sect. 4.1 one can estimate the momentum rate of 
the "J2" outflow from the average maser distance and veloc- 
ity. Using the peak of the SMA 1.3 mm continuum (see Fig. [7] 
lower panel), which, within the measurement error, is coinci- 
dent with the VLA source "A", as the best-guess position of the 
MYSO, and considering only the "J2-SW" water masers, trac- 
ing a collimated flow, the average distance and speed are 0762 
and 38 km s _1 , respectively. Taking the water masers of the 
"J2-NE" group, with more scattered direction of motion, the av- 
erage distance and speed are 171 and 19 km s _1 . In deriving 
these values, we have averaged over a similar number of per- 
sistent maser features, 12 and 13 for the southwestern and 
northeastern cluster, respectively. For the "J2-NE" maser group, 
the separation from the MYSO is about twice as great as that 
of the masers in the "J2-SW" cluster, and the average speed is 
approximately half as much. Then, using Eq. ([T]i for the jet mo- 
mentum rate, the values determined using the southwestern and 
northeastern maser clusters are in good agreement. That lends 



support to our argument that the same jet is driving the motion 
of both "J2" maser lobes. The momentum rate of the "J2" jet is 
P ~ 50(Q/47r) M q yr -1 km s _1 , where Q. is the solid angle of 
the jet. 

In an attempt to search for outflows emerging from the HMC 
in G31.41+0.31, |Cesaroni et al.| ( 12011) have used the SMA to 
map the emission of the CO(2-l) line with an angular resolution 
of «078. The high-velocity (10 km s _1 < |V LS r- V sys | <20 km 
s _1 ) 12 CO(2-l) emission presents a bipolar structure, with red- 
and blue-shifted emission peaks offset to the East and the West 
of the HMC, respectively (see Cesaro nTet al.|20lT[ Fig.l 1). The 
line crossing the two 12 CO(2-l) emission peaks forms a small 
angle (<20°) with the elongation axis (at PA=s70°) of the "J2" 
water maser distribution. The separation on the sky of the two 
emission lobes, «5", is larger than the extension, ~2", of the 
"J2" maser group along its major axis. That suggests that the wa- 
ter masers in "J2" trace the root of a collimated outflow which 
has propagated throughout the HMC accelerating the core gas up 
to 20 km s _1 (along the line-of-sight) at separations of ~10 4 AU 
from the powering MYSO. That the observed bipolar (East- 
West) distribution of the high-velocity 12 CO(2-l) emission could 



mark a col limated outflow has been also discussed by Cesaroni 



etal. 



(2011). Using the CH' 3 CN(12-11) lines, which show a well 



defined Vlsr gradient (at PA = 68°) consistent in orientation and 
amplitude with that seen in I2 CO(2-l) and, in addiction, being 
optically thin, permit a reliable estimate of the whole flow mass, 
these authors derive a value for the outflow momentum rate of P 
= 0.3 M Q yr 1 km s -1 . From the expression quoted above for the 
momentum rate of the jet driving the "J2" water masers, such 
a jet would require only a small solid angle Q. ^0.07 sterad 
(corresponding to a semi-opening angle 9j ~9°), to be powerful 
enough to sustain the molecular flow observed at larger scales. 

The momentum rate derived for both the "Jl" and "J2" water 
maser jets is of the order of a few tenths of M Q yr~' km s _1 . 
Single-dish studies with 10" angular resolution of molecular 
outflows in massive star-forming regions indicate that such val- 
ues of momentum rates are close to the upper limits measured 
for outflows from MYSOs (Lopez-Sepulcre et al. 2009 and 
refs. therein). On the other hand, interferometric observations at 
sub-arcsecond angular resolution have identified a small sample 
of MYSOs (HH 80-81 (parti et al.|[l998); IRAS 16547-4247 



(Rodrigue z et al.| |2008); IRAS 20126+4104 ( |Moscadelli etal. 
201 1 1) which, although being of moderate luminosity (bolo- 
metric luminosity in the range 1 4 — 1 5 L ), eject very power- 
ful, collimated thermal jets and molecular outflows (momentum 
rates of the order of 0.1 M yr -1 km s" 1 ). These results suggest 
that the most compact and youngest (dynamical time <10 4 yr) 
outflows in massive star-forming regions (observable only with 
sub-arcsecond angular resolution) are also among the most pow- 
erful and collimated ones observed. Note that the dynamical time 
scale derived for the "Jl" jet using the water masers is only of 
1.3 10 3 yr, and that of the "J2" jet deduced from the (s patial and 
velocity) separation of the 12 CO(2-l) lobes is 4 10 3 yr (Cesaroni 
|et al |201 1| >. 



4.3. Methanol and OH Maser Dynamics: Jet-driven 
Expansion of the HMC 



As noted in Sect. |3.3| from the comparison of the positions and 
proper motions of the 6.7 GHz methanol (Fig. |4) and 22 GHz 
water (Fig.|2]i masers, one sees that the two clusters of 6.7 GHz 
masers near the "Jl" jet (see Fig. [5] central and right pan- 
els) show consistent velocities with the adjacent water masers, 
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despite the lower speeds. The positions and velocities of the 
6.7 GHz masers near the "Jl" jet suggest that this maser emis- 
sion also traces gas entrained in the jet. To explain the lower ve- 
locities, one could argue that the 6.7 GHz masers would emerge 
from gas at larger separation from the jet axis than the wa- 
ter masers. Concerning the overall distribution of the 6.7 GHz 
masers, Figure |4] (upper panel) evidences that most of them are 
found within a few tenths of arcsec (or a few 1000 AU) from 
the "Jl" water masers (i.e. inside the two dotted rectangles of 
Fig. HJ; in particular, 6.7 GHz masers are also detected at in- 
termediate positions between the "Jl-S" and the "J2-SW" clus- 
ters (see Fig. [T| The proximity to the jets may indicate that the 
6.7 GHz maser excitation could be favoured by the outflow ac- 
tivity. One can speculate that the mid-infrared photons, required 
for pumping the 6.7 GHz masers (Cragg et al. 2002). could more 
easily penetrate into the dense HMC gas in proximity of the out- 
flows, taking advantage of the cavities induced by the outflow- 
ing gas. The proposed association of the 6.7 GHz masers with 
outflows, based on the similarity of the 3-D velocities of the 



methanol and water masers, agrees with the results by De Buizer 
(2003| and |De Buizer et a l. (2009), which find that in most cases 
the methanol masers are elongated parallel to jets traced by the 
H2 and SiO emissions, suggesting a physical association of the 
6.7 GHz masers with the corresponding outflows. 

Looking at the overall distribution of 6.7 GHz maser ve- 
locities, Fig. |4] (upper panel) shows that most of the 6.7 GHz 
maser proper motions are diverging from the HMC center, with 
minor deviations from the radial direction. We have assumed 
the HMC center to be the midpoint between the VLA sources 
"A" and "B". Altogether, the pattern of 6.7 GHz proper mo- 
tions marks the expansion (at an average speed of »10 km s _1 ) 
of the HMC gas close to the "Jl" and "J2" jets. Observing the 
G31. 41+0.31 HMC with SMA at sub-arcsecond angular resolu- 
tion, |Girart et al.| (|2009 ) detected an inverse P-Cygni profile in 
the spectrum of the C S (7-6) line, suggesting infall at veloc- 
ity of a few km s _1 of the core envelope at radii of «10 4 AU. 
Although the 6.7 GHz masers trace comparable radial distances, 
in the range 10 3 — 10 4 AU, the pattern of their 3-D velocities 
presents very limited clues to infall. It may be significant that the 
only inward motions of the 6.7 GHz masers are observed for the 
two features with the largest distance from the water maser jets 
(inside the dotted ellipse in Fig. [4] upper panel). If the 6.7 GHz 
masers are excited mostly in proximity of the "Jl" and "J2" jets, 
their average, outward motion might not be representative of the 
motion of the HMC bulk gas, since the gas kinematics at large 
separation from the jets is not adequately sampled. 

Figure[3]illustrates that the OH 1 .6 GHz masers have an East- 
West elongated distribution with an enough regular increase of 
maser Vlsr with position varying from West to East. The Vlsr 
gradient traced by the OH masers has a similar orientation as the 
C H 3 CN Vlsr gradient a cross the G3 1.41 +0.31 HMC revealed 
by |Cesaroni et al.| ( |201 1\ . Figure |6]presents the OH maser Vlsr 
versus the corresponding positions projected along the direc- 
tion of the CH3CN gradient (PA = 68°). The maser velocities 
are well correlated with the positions (correlation coefficient = 
0.8) with a gradient of 4.7±1.6 km s _1 arcsec -1 . Inspecting the 
position-velocity plot of the 12 CO(2-l) line, observed with the 
SMA by Cesaroni et al.| ( |201 1) , one notes that a similar (both 
in orientation and in amplitude) Vlsr gradient is traced by the 
high-velocity wings of the 12 CO(2-l) emission (offset in pro- 
jected position by 4"-5" in correspondence of the maximum 
Vlsr separatio n of 20-30 km s _1 , see Fig. [6J. Since, as dis- 



"J2" water maser jet, we postulate that the observed OH 1 .6 GHz 
masers originate from the same outflow. The similarity of the 
Vlsr gradient detected in the high-velocity 12 CO(2-l) emission 
wings and in the OH masers can indicate that the OH masers are 
excited relatively close to the outflow axis (with respect to the 
transverse size of the outflow), where the HMC gas is pushed by 
the jet to relatively higher velocities. That could explain why in 
Fig.|6]the OH maser Vlsr are offset from those of the CH3CN, 
which should trace the velocity of the HMC bulk gas. 

It is interesting to note that, whereas the "Jl" jet appears to 
be responsible for the excitation and the motion of most of the 
6.7 GHz masers (as discussed above), the OH 1.6 GHz masers 
seem to be mainly associated with the "J2" jet. These two maser 
types are both radiatively excited and models predict similar 
physical excitation conditions (Cragg et al. 2002 ), except for rel- 
atively high gas temperature (=s!50 K), and densities (10 7 cm -3 



< «h 2 ^ 10 cm ), for which OH 1.6 GHz masers are pre- 



dicted to thermalize. From CH3CN observations, Beltran et al 



(2005| derive a temperature and a density for the HMC gas of 
^100 K and ~10 7 cm -3 , which is close to the predictions of the 
models for OH 1.6 GHz maser quenching. In Sect. 



4. 1 on the 



cussed in Sect. 4.2 the CO(2-l) line is likely excited in the 



collimated, East-West oriented, molecular outflow driven by the 



basis of the size and the shape of the water maser distribution, 
we have postulated that the "Jl" jet has not yet emerged from the 
MYSO natal cocoon, which could be too dense and warm for ef- 
ficient pumping of the OH 1.6 GHz masers. On the other hand, 
the more extended and evolved "J2" jet could have blown the 
densest portions of the HMC gas away from its axis, reducing 
the average gas density to «h 2 < 10 7 cirT 3 , suitable for exciting 
the OH 1.6 GHz masers. A reduced gas density, however, can- 
not explain the relatively few 6.7 GHz masers observed towards 
the "J2" jet, since this maser emission should be strong even at 
densities $10 7 crrT 3 . One possibility is that the environments of 
the "Jl" and "J2" jets present chemical differences affecting the 
abundance rates between CH3OH and OH. 

5. Discussion 

5. 1 . "J1" and "J2" Jets and the V L sr Gradient across the 
HMC 

One of the most notable findings of previous interferometric 
(PdBI, SMA) observations towards G3 1.41 +0.31 is the discov- 
ery of a well defined Vlsr gradient traced in several transitions 
of high-density tracers (remarkably CH3CN, see Fig. [7} and also 
in the 12 CO(2-l) line, a typical outflow tracer. The interpreta- 
tion of this velocity gradient, in terms either of a rotating toroid 
or a collimated outflow, is still under debate. In the following, 
we discuss the relevance of the VLBI observations of the water 
masers to shed light on this issue. 

5.1 .1 . The Origin of the "J1 " Jet 

One of the main difficulty to interpret the Vlsr gradient in 
terms of (pseudo-)Keplerian rotation is the absence of a jet ori- 
ented perpendicular to the Vlsr gradient, as predicted by star- 
formation models. We have detected however the "Jl" water 
maser distribution, which traces a collimated jet traveling along 
a direction forming an angle of «60° with the gradient orienta- 
tion. According to the interpretation by Cesaroni et al. ( 201 1) , 
the G3 1.41 +0.31 HMC is a massive (several 100 M G ), large ro- 
tating toroid feeding a cluster of MYSOs embedded at the cen- 
ter of the core. One could expect that the accrection disks of 
these MYSOs are aligned along directions at close angle with 
the Vlsr gradient orientation. Then it is plausible that the "Jl" 



12 



L. Moscadelli et al.: A Double- Jet System in the G3 1.4 1+0.31 Hot Molecular Core. 



CO 

E 



120 



110 



100 



01 

in 



> 



90 - 



80 - 



12 



CO, CH 3 CN, OH 




NE 



-2 

Offset (arcsec) 



SW 



Fig. 6. Position-velocity plots of the 12 CO(2-l) (grey scale) and CH 3 CN(12-11) K=4 line emission obtained with the SMA by Cesaroni et al. 
(201 1). The offset is measured along the direction with PA = 68° passing through the HMC (see Fig. [3j and is measured from the phase center of 
the SMA observations (o-(J2000)=18 h 47 m 34 ! :315, (5(J2000)=-01° 12'45'.'9), assumed positive toward NE. Contour levels increase from 0.2 in steps 
of 0.2 Jy beanT 1 for CH 3 CN and from 0.24 in steps of 0.24 Jy beairT 1 for 12 CO. The cross in the bottom left denotes the angular and spectral 
resolutions of the SMA data. The solid dots represent the 1665 MHz and 1667 MHz OH maser features detected in our VLBA observations. The 
black dashed line gives the best linear fit of the OH maser Vlsr versus corresponding positions. . 



jet, forming an angle of «30° with the direction perpendicular to 
the Vlsr gradient, is associated with one of these MYSO. |Osorio| 
et al.| ( |2009] > have recently reproduced the dust spectrum and the 



ammonia line emission of the G31.41+0.31 HMC with a model 
of an infalling envelope onto a massive star undergoing an in- 
tense accretion phase, deriving a central star mass of 20-25 M Q 
and an age of 3^4 10 4 yr. Note that the modeled protostar age 
is more than one order of magnitude larger than the dynamical 
time derived for the "Jl" jet of 1.3 10 3 yr. That would lead us to 
rule out the hypothesis that the "Jl" jet is associated with the ac- 
creting protostar, unless our estimate of the "Jl" dynamical time 
is not severely underestimated. 

If the "Jl" jet is precessing, that would explain both the 
underestimation of the jet time scale and why the "Jl" jet 
is not oriented perpendicular to the Vlsr gradient. In the lit- 
erature, it exists at least one, well-studied case, the MYSO 
IRAS 20126+4104, which emits a jet precessing at an angle 
as large as 45° (Shephe rd et al.p OOO). The observation of the 
two VLA compact sources "A" and "B" at the center of the 
G31. 41+0.31 HMC might be hinting at a massive binary sys- 
tem, and the tidal interaction between the disk and a compan- 
ion star in a non-coplanar orbit is one possible mechanism for 
jet precession (Terquem et al. 1999 Bate et al. 2000} . We do 
not find any hint of a jet along "Jl" in the SMA CO maps by 
Cesaroni et al. ( 201 1) . The present data towards G31.41+0.31 



do not allow to discuss the case for precession in a more quan- 
titative way. Future ALMA observations in thermal continuum 
and line emissions, achieving an angular resolution ^ O'.'l, could 
permit to identify the circumstellar disk around each component 
of the (putative) binary system, to verify if the disk planes are 



misaligned with respect to the binary orbit, and to determine the 
binary separation. 



5.1.2. "J2" Jet Interaction with the HMC 

The "J2" water maser distribution (see Figs.[T]and[2} is elongated 
along a direction (PA ^70°) in good agreement with the orien- 
tation (PA = 68°) of the V LS r gradient dete cted in G3 1.41 +0.3 1 
with the SMA in different lines of CH 3 CN ( |Cesaroni etal.|2011) . 
When discussing the outflow interpretation for the Vlsr gra- 
dient, Cesaroni et al. ( 201 1 1 derive a momentum rate of P — 
0.3 M Q yr" 1 km s" 1 from the optically thin CH' 3 CN(12-1 1) line. 
This value is consistent with our estimate of P for the "J2" wa- 
ter maser jet (see Sect. |4.2[> for a beaming angle Q »0.07 sterad. 



We thus favour the interpretation of the CH3CN Vlsr gradient in 
terms of a compact and collimated outflow. 



Cesaroni et al. (201 \) point out two major difficulties for the 



outflow interpretation of the Vlsr gradient. The first one is asso- 
ciated with the too large flow parameters (in particular the mo- 
mentum rate) if compared with the results of single-dish surveys 
of massive outflows. As already discussed in Sect. |4.2| high- 
angular resolution observations of a few, well studied MYSOs 
suggest that, analogously to what observed in low-mass pro- 
tostars, the youngest and most compact outflows of massive 
(proto)stars are also intrinsically more powerful. Thus, the re- 
sults of single-dish surveys, biased towards more evolved out- 
flows, could systematically underestimate the flow parameters 
at the earliest evolutionary phases (dynamical time < 10 4 yr). 
The second concern of Cesaroni et al. ( 201 1 1 is that the flow dy- 
namical time (4 10 3 yr) would be too short with respect to the 
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time needed to form typical hot core species - such as methyl 
cyanide - according to some theoretical models (Charnley et al. 
1992). The formation route of these models requires that the gas 



spends at least several 10 3 yr at temperatures above 100 K, since 
the process is powered by the evaporation of H2O and other 
non-refractory species at 100 K, followed by a chain of high- 
temperature gas-phase reactions to transform simple carbon- 
bearing molecules into second-generation complex organics. 
However, these high-temperature models are presently chal- 
lenged by the grain-surface models, which predict that the or- 
ganics may be formed abundantly on the grain surfaces through 
mild p hotochemistry at 20^-0 K dGarrod & Herbst 2006}|Garrod 



2008). The grain-surface pathway takes several 10** yr at 20- 



40 K, and such timescales and temperatures are attainable both 
in the envelope before the material approaches the star and inside 
the disk (should it be present). Once formed on the grains, these 
first-generation organics can be returned to the gas phase by ther- 
mal evaporation or liberated by non-thermal processes such as 
shocks in the outflow walls. The whole formation route requires 
that the material spends only a few 100 yr at high-temperatures, 
an order of magnitude less than the dynamical time scale of the 
outflow supposed to be at the origin of the observed Vlsr gradi- 
ent in the G31.41+0.31 HMC. 

Although the observed water maser kinematics is consis- 
tent with the outflow scenario, the maser data themselves are 
not clear-cut enough to rule out alternative interpretations. In 
Sect. |4.2| we noted that the water maser velocities to the East of 
the VLA sources "A" and "B" are not collimated along the axis 
of the "J2" maser distribution, which might cast some doubts on 
these masers to belong to a collimated outflow, and on the ef- 
ficacy of the "J2" flow to produce a well-defined Vlsr gradient 
across the HMC. We postpone the discussion of this problem to 
Sect. |5.2[ where we show how also the anomalous velocities of 
these features can be explained in the contest of the jet interpre- 
tation for the "J2" masers. 

5.2. Interaction of the HMC with its Environment 

In the upper panel of Fig. [JJ we compare the velocity map of 
the HMC obtained from the CH3CN molecule with the free-free 



continuum emission imaged at 1.3 cm by |Cesaroni et al. ( 1994b 
2010). Two continuum maps are shown: one obtained by merg- 
ing the C and D configurations of the VLA, the other with the B 
configuration. Two peaks of emission are seen in the former im- 
age, with the one to the SE being located very close to the HMC. 
In the B-array map only the peak to the NE is still visible: this 
is very likely tracing the position of the 06 star (2 x 10 3 L e ; 
Cesaroni et al. (2010)) responsible for the extended free-free 
emission seen in Fig. [7j upper panel. The UV photons escap- 
ing from the surroundings of the star could also ionize the sur- 
face of the HMC, thus causing the second, more shallow peak 
of free-free continuum observed with the C+D arrays. This sce- 
nario strongly suggest that the star ionizing the HII region and 
the HMC are close in space and not only in projection on the 
plane of the sky. Further support to this hypothesis is lent by 
the shape of the millimeter continuum emission, shown in the 
bottom panel of Fig. [7J Clearly, the contour map looks slightly 
elongated in the NE-SW direction, with the emission extending 
towards NE and decreasing more sharply to the SW. Such an 
asymmetry can be explained with heating of dust on the side of 
the HMC facing the luminous star ionizing the HII region. 

We believe that the proximity between the ionizing star and 
the HMC, besides explaining the morphology of both the free- 
free and thermal dust continuum emission, can also justify the 



directions of proper motions in the "J2-NE" H2O maser group. 
If the "J2" masers are tracing a bipolar outflow, the NE lobe 
of this can be strongly affected by the interaction with the stel- 
lar photons, possibly enhancing the level of turbulence of the 
outflowing molecular gas through heating and photoevaporation 
( Henney et al.|2009) . Albeit impossible to quantify on the basis 
of the available data, the effect of these phenomena on the ve- 
locity field of the outflowing gas is likely significant and could 
explain the anomalous maser proper motions observed. 

It is also interesting to note that a significant enhancement 
of turbulence could also explain why no CH3OH maser and 
only few OH masers are detected on the NE side of the HMC. 
These types of masers have measured diameters ^ 10 mas, much 
greater than those of H2O masers (~ 0.5mas), which suggests 
that also the typical amplification paths are larger. Therefore, 
OH and CH3OH masers are also those most affected by an al- 
teration of the gas velocity field, consistent with the observed 
lack of emission. 

In conclusion, we believe that scattered proper motions of 
the "J2-NE" cluster of water masers cannot be used to discard the 
hypothesis that such masers are tracing the NE lobe of a bipolar 
outflow oriented NE-SW. 



6. Conclusions 

This work reports on multi-epoch VLBI observations of the H2O 
22 GHz and CH3OH 6.7 GHz masers, and single-epoch VLBA 
of the OH 1.6 GHz masers towards the HMC in G3 1.4 1+0.31. 
Distributed over the HMC, we detect 173 H 2 0, 85 CH 3 OH, and 
69 OH maser features; proper motions are measured for 70 H2O 
and 40 CH3OH, persistent features. The magnitude of the wa- 
ter maser, relative proper motions ranges from =sl0 km s _1 to 
« 120 km s _1 , with a mean value of 34.0 km s , and a mean 
error of 7.0 km s . The magnitude of the methanol maser, rela- 
tive proper motions ranges from ^5 km s _1 to a;30 km s _1 , with 
a mean value of 13.0 km s _1 , and a mean error of 4.0 km s . 

Water masers present a symmetric spatial distribution with 
respect to the HMC center, where two nearby (072 apart), 
compact, VLA sources (named "A" and "B") are detected by 
2010p . A first group of water masers, named 



Cesaroni et al. 



"Jl", consists of two bow-shape distributions, offset to the North 
and the South of the VLA sources by ^075, with the bow tip 
pointing to the North (South) for the northern (southern) cluster. 
The "Jl" distribution is well fit with an elliptical profile (ma- 
jor and minor axis of l'.'44±0'.'l and 0'.'24±0'.'02, respectively, 
with the major axis PA (East from North) = 8°±1°), and the 
maser proper motions mainly diverge from the ellipse center, 
with average speed of 36 km s . These findings strongly sug- 
gest that the "Jl" water maser group traces the heads of a young 
(dynamical time of 1.3 10 3 yr), powerful (momentum rate of 
^0.2 M Q yr~' km s~'), collimated (semi-opening angle ^10°) 
jet emerging from a MYSO located close (within «0'.'15) to the 
VLA source "B". 

Most of the water features not belonging to "Jl" present 
a NE-SW oriented, S -shape distribution, which consists of an 
NE-SW extended (^2"in size, PA«70°) strip of features, and 
two smaller (=s075 in size), SE-NW oriented lines of masers 
at the ends of the NE-SW strip. We denote this S-shape, wa- 
ter maser group with "J2". The few measured proper motions 
for water masers belonging to the NE-SW elongated strip indi- 
cate a motion receding from the two VLA continuum sources 
(located approximately at the center of the strip), along direc- 
tions at close angle with the strip orientation. The water features 
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by steps of 0.2 Jy beam . The white dots denote the two VLA compact sources detected by ( 
vectors give the absolute positions and relative proper motions (with respect to the "center of motion", see Sect. |3. 1 1 > of the 22 GHz water masers. 
Triangle area is proportional to the logarithm of the maser intensity. The amplitude scale for proper motions is indicated by the black arrow in the 
bottom left corner of the panel. The S-shape dashed polygon encompasses the "J2" water maser distribution and the dotted circle surrounds the 
isolated cluster of water masers belonging to neither the "J2" nor the "Jet 1" group. 
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distributed at the southwestern "J2" end move all to the SW, par- 
allel to the strip orientation. Instead, at the northeastern end of 
the strip, the velocities of the water maser features point in dif- 
ferent directions. The elongated distribution of the "J2" group 
and the direction of motion, approximately parallel to the di- 
rection of elongation, of most "J2" water masers suggests that 
the masers are tracing another collimated outflow, emitted from 
a MYSO placed near the VLA source "A". The more scattered 
distribution of proper motions of the water masers at the "J2" 
northeastern end could reflect turbulence induced by the heat- 
ing and photoevaporation of the northeastern side of the HMC, 
which is facing the O-type star responsible for the H II region 
observed 5" to NE of the HMC. 

The orientation (PA w 70°) of the "J 2" jet agrees well with 
that (PA = 68°) of the V LSR gradient across the HMC revealed in 
several transitions of CH3CN and also in the 12 CO(2-l) line by 
|Cesaroni et al.| ( [201 Besides, the "J2" jet is powerful enough 
to explain the large mo mentum rate, . 3 M„ yr" 1 km s , esti- 



CN by 



Cesaroni et al. 



(201 1 ) under the hypoth- 



mated from CH^ 3 
esis that the Vlsr gradient represents a collimated outflow. These 
two facts lead us to favour the interpretation of the CH3CN 
Vlsr gradient in terms of a compact and collimated outflow. The 
proper motions of the CH3OH 6.7 GHz masers, mostly diverging 
from the HMC center, also witness the expansion of the HMC 
gas driven by the "Jl" and "J2" jets. 

As already indicated by previous results for other MYSOs, 
this work on the G3 1.4 1+0.31 HMC confirms that VLBI maser, 
3-D kinematics can be relevant for the study of massive star- 
formation. This technique is presently the only one which en- 
ables to resolve the typical sizes (several 100 AU) of massive 
stellar clusters, and to study in 3-D the process of mass accre- 
tion and ejection around a single MYSO. 
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L. Moscadelli et al.: A Double- Jet System in the G3 1.4 1+0.31 Hot Molecular Core. 
Table 2. 1665 MHz OH Maser Parameters. 



Feature 


Pol. 


Ipeak 




Ax 


Ay 


Number 




(Jy beam -1 ) 


(km s" 1 ) 


(mas) 


(mas) 


1 


LCP 


4.85 


98.0 


0.00±0.00 


0.00+0.00 


2 


LCP 


0.86 


111.4 


1906.92±0.35 


160.94+0.47 


3 


LCP 


0.57 


103.7 


695.99±0.70 


380.96+1.05 


4 


LCP 


0.43 


93.7 


552.01±0.77 


-420.45+0.94 


5 


LCP 


0.41 


94.8 


709.72±0.70 


-560.67+0.92 


6 


LCP 


0.23 


102.9 


692.67±0.99 


330.48+1.41 


7 


LCP 


0.16 


96.9 


241.37±1.34 


-42.28+1.94 


8 


LCP 


0.15 


102.2 


-317.18±1.50 


-80.13+1.85 


9 


LCP 


0.14 


97.9 


-692.10±1.84 


-53.60+2.05 


10° 


LCP 


0.11 


109.3 


-1757.31±2.26 


-842.74+2.71 


1 


RCP 


1.49 


94.8 


711.06±0.26 


-560.97+0.33 


2 


RCP 


0.82 


97.7 


-690.91+0.47 


-46.86+0.56 


3 


RCP 


0.76 


96.0 


549.14+0.40 


-424.47+0.49 


4 


RCP 


0.73 


97.6 


-769.73+0.52 


-67.21+0.55 


5 


RCP 


0.41 


103.2 


696.07+0.82 


374.85+1.12 


6 


RCP 


0.32 


102.8 


695.33+1.13 


312.30+1.24 


7 


RCP 


0.24 


97.6 


-2.33+1.81 


-4.39+2.11 



Notes. Column 1 gives the feature label number; column 2 indicates the circular polarization of the maser emission; columns 3 and 4 report the 
intensity of the strongest spot and the intensity-weighted LSR velocity, respectively; columns 5 and 6 list the position offsets (with the associated 
errors) along the R.A. and Dec. axes, relative to the feature #1 of LCP. 

{a) This single feature is significantly offset to West with respect to the other OH maser detections, and it is not plotted in Figure|3] 
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L. Moscadelli et al.: A Double- Jet System in the G3 1.4 1+0.31 Hot Molecular Core. 
Table 3. 1667 MHz OH Maser Parameters. 



Feature 


Pol. 


Ipeak 


Vlsr 


A x 


Ay 


Number 




(Jy beam -1 ) 


(km s ') 


(mas) 


(mas) 


1 


LCP 


6.20 


97.8 


0.00±0.00 


0.00±0.00 


2 


LCP 


3.42 


97.9 


-777.01±0.19 


-69.90±0.19 


3 


LCP 


2.30 


103.6 


702.06±0.21 


371.33±0.23 


4 


LCP 


1.65 


96.2 


-713.49±0.29 


-128.56±0.25 


5 


LCP 


1.61 


103.8 


731.84±0.25 


411.39±0.29 


6 


LCP 


1.08 


95.3 


-743.87±0.38 


-132.20+0.32 


7 


LCP 


0.99 


96.9 


-637.47±0.38 


-122.11+0.40 


8 


LCP 


0.94 


111.2 


1907.56±0.46 


164.47+0.49 


9 


LCP 


0.72 


108.7 


1778.70±0.62 


84.73+0.57 


10 


LCP 


0.71 


89.5 


-372.31±0.48 


-540.98+0.49 


11 


LCP 


0.53 


97.3 


-669.36±0.56 


-43.81+0.62 


12 


LCP 


0.41 


97.0 


-810.71±0.77 


-54.05+0.71 


13 


LCP 


0.37 


96.6 


239.25±0.85 


-39.08+1.04 


14 


LCP 


0.27 


92.1 


-961.22±1.22 


-68.13+1.30 


15 


LCP 


0.26 


103.1 


545.39±1.12 


404.28+1.18 


16 


LCP 


0.25 


100.3 


-446.38±1.16 


-200.60+1.26 


i / 


LLP 


U.24 


101.9 


/2U. 38+1.4/ 


20y. / 0+1. 43 


18 


LCP 


0.23 


99.0 


-70.31±1.10 


-29.90+1.16 


19 


LCP 


0.22 


92.9 


-1057.53±1.71 


-22.39+1.63 


20 


LCP 


0.17 


102.0 


561.21+1.73 


141.02+1.91 


21 


LCP 


0.17 


96.0 


-316.07±1.63 


-332.44+1.97 


22 


LCP 


0.15 


94.1 


545.39±1.42 


-426.44+1.68 


23 


LCP 


0.15 


93.2 


-909.66±1.86 


-87.36+1.59 


24 


LCP 


0.12 


100.6 


-206. 02± 1.84 


-768.08+2.36 


1 


RCP 


6.13 


97.6 


-773.llzt0.17 


-69.95+0.16 


2 


RCP 


2.92 


103.5 


723.50±0.28 


399.02+0.31 


3 


RCP 


2.25 


103.2 


700.38±0.20 


366.62+0.22 


4 


RCP 


1.28 


96.3 


-715.28±0.42 


-129.39+0.35 


5 


RCP 


1.09 


97.6 


-698.57±0.62 


-50.59+0.54 


6 


RCP 


0.91 


97.3 


-6.42±0.53 


-2.40+0.55 


7 


RCP 


0.68 


97.6 


-740.27±1.11 


-615.34+0.87 


8 


RCP 


0.63 


96.8 


-637.96±0.61 


-119.64+0.60 


9 


RCP 


0.61 


97.8 


-740.02+0.78 


-62.81+0.82 


10 


RCP 


0.55 


92.0 


-961.68±0.72 


-65.10+0.66 


11 


RCP 


0.54 


97.1 


-671.03±0.59 


-45.12+0.60 


12 


RCP 


0.46 


95.4 


-747.35±0.82 


-131.52+0.73 


13 


RCP 


0.44 


95.7 


240.08±0.69 


-38.21+0.84 


14 


RCP 


0.41 


96.8 


-809.27±0.94 


-52.88+0.89 


15 


RCP 


0.38 


102.8 


694.81±0.94 


311.94+0.97 


16 


RCP 


0.37 


98.9 


-70.41±0.91 


-29.06+0.94 


17 


RCP 


0.35 


111.1 


1907.42±0.75 


169.23+0.92 


18 


RCP 


0.34 


89.9 


-372.86±0.72 


-542.73+0.75 


19 


RCP 


0.34 


96.6 


-671.38ztl.39 


438.76+1.16 


20 


RCP 


0.32 


102.4 


-1196.87±1.14 


-50.74+1.33 


21 


RCP 


0.31 


95.3 


-321.55±0.91 


Q 1 1 AC\ i 1 f\A 

-331.49+1.04 


22 


RCP 


0.27 


93.3 


-917.20±1.19 


-89.33+1.20 


23 


RCP 


0.23 


96.6 


-675.39±1.54 


-116.40+1.32 


24 


RCP 


0.23 


98.9 


644.40±1.10 


-272.34+1.34 


25 


RCP 


0.17 


95.4 


547.94±1.78 


-424.08+2.35 


26 


RCP 


0.15 


102.9 


546.51 ±2.07 


412.25+2.21 


27 


RCP 


0.15 


91.1 


-901.27±1.60 


-75.20+1.66 


28 


RCP 


0.13 


89.8 


-706.07±1.90 


-182.16+2.37 



Notes. Column 1 gives the feature label number; column 2 indicates the circular polarization of the maser emission; columns 3 and 4 report the 
intensity of the strongest spot and the intensity-weighted LSR velocity, respectively; columns 5 and 6 list the position offsets (with the associated 
errors) along the R.A. and Dec. axes, relative to the feature #1 of LCP. 
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L. Moscadelli et al.: A Double- Jet System in the G3 1.4 1+0.31 Hot Molecular Core. 
Table 4. 22 GHz H,0 Maser Parameters 



Feature 


Epochs of 


Ipeak 


Vlsr 


A x 


Ay 


v x 


V y 


Number 


Detection 


(Jy beam -1 ) 


/l - 1 \ 

(km s 1 ) 


(mas) 


(mas) 


(km s 1 ) 


/l -1 \ 

(km s 1 ) 





1,2,3,4 




98.4 


173.65±0.06 


767.21+0.06 


0.0+0.0 


0.0+0.0 


1 


1,2,3,4 


108.04 


98.0 


0.00±0.00 


0.00+0.00 


0.6+5.7 


-30.7+5.9 


2 


1,2,3,4 


30.13 


97.5 


760.41±0.08 


211.81+0.08 


-13.3+5.7 


1.3+5.8 


3 


1,2,3,4 


21.02 


97.2 


-625.59±0.08 


390.26+0.08 


-12.4+5.4 


-6.4+5.5 


4 


1,2,3,4 


20.30 


100.9 


-16.62±0.07 


-5.71+0.08 


6.4+5.5 


-36.0+5.7 


5 


1,2,3,4 


12.75 


95.2 


-514.10+0.08 


280.57+0.08 


-24.9+5.6 


-29.1+5.7 


6 


1 


11.49 


97.8 


759.52+0.07 


211.05+0.08 






7 


1,2,3 


8.99 


94.8 


123.56+0.08 


476.60+0.08 


36.5+9.1 


1.9+9.5 


8 


1,2,3,4 


6.71 


100.1 


827.44+0.07 


1418.08+0.08 


16.0+5.5 


0.1+5.6 


9 


1,2,3,4 


6.45 


93.6 


-514.71+0.08 


247.81+0.08 


-21.8+5.6 


-7.2+5.7 


10 


1,2,3,4 


6.41 


97.0 


772.36+0.08 


245.23+0.08 


-2.1+5.5 


26.3+5.7 


11 


1,2,3,4 


6.37 


99.5 


15.45+0.08 


14.27+0.08 


-3.1+5.5 


-18.6+5.6 


12 


1,2,3,4 


5.83 


102.1 


860.39+0.08 


1371.82+0.08 


3.7+5.5 


19.4+5.6 


13 


1,2,3,4 


5.51 


98.3 


-2.99+0.07 


0.20+0.08 


65.6+5.5 


-21.4+5.8 


14 


1,2,3,4 


4.95 


96.0 


-513.87+0.07 


282.77+0.08 


-20.5+5.4 


-12.4+5.5 


15 


3 


4.71 


96.0 


-12.21+0.07 


-3.91+0.08 






16 


1,3 


4.67 


97.5 


761.64+0.07 


213.02+0.08 






17 


1,2,3,4 


4.66 


100.6 


596.73+0.08 


1444.51+0.08 


-4.1+5.7 


15.0+5.9 


18 


1,2,3,4 


4.28 


99.8 


823.47+0.07 


1424.87+0.08 


5.6+5.5 


-0.7+5.6 


19 


3 


4.27 


87.4 


106.77+0.07 


458.18+0.08 






20 


1,2,3,4 


4.20 


113.2 


-106.13+0.08 


432.19+0.08 


28.9+5.8 


-60.0+6.0 


21 


1,2,3,4 


4.19 


98.5 


10.63+0.08 


5.19+0.08 


40.5+5.5 


-28.3+5.7 


22 


1,2,3,4 


3.89 


97.2 


771.83+0.07 


239.24+0.08 


1.3+5.4 


26.6+5.7 


23 


1,2,3,4 


3.77 


99.3 


884.14+0.07 


1345.64+0.08 


14.7+5.5 


18.7+5.6 


24 


1,2,3,4 


3.72 


97.6 


771.01+0.08 


234.36+0.08 


-2.5+5.5 


23.6+5.7 


25 


1,2,3 


3.13 


118.4 


41.31+0.08 


1187.73+0.08 


-9.7+9.0 


82.7+9.2 


26 


1,2,3,4 


3.12 


99.3 


881.46+0.08 


1326.57+0.08 


10.0+5.5 


11.5+5.8 


27 


3 


2.43 


93.1 


-12.38+0.07 


-3.67+0.08 






28 


1,2,3,4 


2.30 


100.2 


207.47+0.08 


1228.88+0.08 


10.3+5.5 


-4.4+5.8 


29 


2,3,4 


2.28 


97.2 


770.52+0.08 


232.86+0.08 


15.8+8.4 


66.7+9.2 


30 


1,2,3,4 


2.28 


89.5 


111.09+0.08 


461.95+0.08 


32.9+5.6 


-50.2+5.8 


31 


2 


1.89 


99.5 


-11.96+0.08 


-3.17+0.08 






32 


1,2,3,4 


1.81 


107.8 


-401.49+0.08 


87.41+0.08 


3.8+6.5 


-79.9+7.1 


33 


1,2,3,4 


1.74 


100.7 


207.63+0.08 


1235.46+0.08 


20.8+5.6 


10.6+5.9 


34 


1,2,3,4 


1.67 


101.2 


-54.19+0.08 


1052.12+0.08 


-17.0+5.4 


26.9+5.7 


35 


2,3,4 


1.62 


114.4 


35.79+0.09 


1184.51+0.11 


-46.6+10.2 


80.2+11.3 


36 


4 


1.55 


87.3 


107.96+0.07 


458.62+0.07 






37 


1,2,3,4 


1.53 


95.8 


-128.10+0.08 


203.66+0.08 


-24.2+5.6 


-4.2+5.8 


38 


1 


1.46 


95.9 


-12.34+0.08 


-3.66+0.08 






39 


1 


1.46 


98.4 


0.74+0.07 


-0.23+0.08 






40 


1 


1.43 


98.7 


-11.99+0.08 


-3.22+0.08 






41 


2,3,4 


1.42 


118.5 


24.08+0.08 


1048.78+0.08 


-12.1+8.7 


45.0+8.9 


42 


2 


1.40 


96.7 


-12.26+0.08 


-3.80+0.08 






43 


1,2,3 


1.36 


96.9 


148.53+0.07 


1412.51+0.08 


8.8+10.4 


1.1+11.7 


44 


1,2,3,4 


1.31 


99.6 


-53.87+0.08 


1057.11+0.08 


-16.3+5.7 


16.6+6.1 


45 


4 


1.30 


96.3 


-12.19+0.08 


-4.02+0.08 






46 


1,2,3,4 


1.29 


94.1 


-8.91+0.08 


1202.27+0.08 


-7.6+5.6 


11.5+5.8 


47 


1,2,3,4 


1.21 


99.3 


970.34+0.07 


1068.50+0.08 


8.4+6.1 


29.7+6.7 


48 


1,2,3,4 


1.19 


95.1 


-514.60+0.08 


275.71+0.08 


-33.8+5.5 


-56.3+5.8 


49 


1 


0.99 


96.6 


15.87+0.08 


8.46+0.08 






50 


1,2,3,4 


0.92 


95.2 


135.06+0.08 


516.27+0.08 


36.6+5.7 


4.6+6.2 


51 


4 


0.77 


94.9 


970.81+0.07 


1069.10+0.08 






52 


3 


0.75 


94.9 


-514.38+0.08 


277.10+0.08 






53 


1,2,3,4 


0.70 


101.8 


10.05+0.08 


1206.54+0.08 


-30.2+5.9 


20.6+6.4 


54 


1,2,3,4 


0.70 


96.7 


103.78+0.08 


472.25+0.08 


44.6+5.9 


14.6+6.1 


55 


1,2,3,4 


0.70 


101.4 


615.77+0.08 


1450.81+0.08 


-8.6+5.7 


17.2+6.0 


56 


1,2,3,4 


0.69 


87.3 


110.86+0.08 


461.48+0.09 


43.9+5.8 


-24.6+6.5 


57 


4 


0.65 


99.7 


-17.09+0.08 


-7.96+0.08 






CO 

58 


1 -"i Q A 

1,2,3,4 


A £i 

0.63 


a*c a 

96.9 


1 cn AC i A AO 

159.05+0.08 


1 1 AC AH i A A A 

1395.47+0.09 


2.5+6.1 


62.3+6.7 


59 


4 


0.61 


92.5 


-12.36+0.08 


-3.79+0.08 






60 


1,2,3,4 


0.59 


96.6 


82.75+0.08 


463.46+0.09 


-3.7+6.1 


-9.3+6.7 


61 


1 


0.56 


101.1 


-405.27+0.08 


89.91+0.08 






62 


2,3,4 


0.55 


103.5 


-406.59+0.08 


90.46+0.09 


-90.0+9.1 


-75.7+9.7 


63 


1,2,3,4 


0.54 


99.6 


878.77+0.08 


1316.82+0.08 


2.9+6.2 


5.7+7.4 


64 


4 


0.54 


98.8 


200.87+0.08 


1309.50+0.08 






65 


1,2,3,4 


0.52 


103.6 


-41.44+0.08 


1125.30+0.08 


-21.4+5.8 


15.6+6.2 



L. Moscadelli et al.: A Double- Jet System in the G3 1.4 1+0.31 Hot Molecular Core. 

Table 4. continued. 



Feature 


Epochs of 


I k 


Vlsr 


A x 


A y 


v x 


V y 


Number 


Detection 


(Jy beam ) 


(km s~') 


(mas) 


(mas) 


(km s~') 


(km s" 1 ) 


66 


3 


0.51 


98.8 


200 78+0 08 


1309 34+0 08 

Uuy . .. ' — !— V./ .V./ y.f 






67 


2 


0.50 


99.9 


-404 50+0 09 


sg 78+0 09 






68 


4 


0.46 


99.3 


15 55+0 08 


8 47+0 08 






69 


4 


0.43 


98.5 


202.48+0.08 


1300 13+0 09 






70 


1,2,3,4 


0.42 


94.7 


103 76+0 08 


1349 41+0 09 


-19.2+5.9 


19.6+6.6 


71 


1,2 


0.42 


100.5 


590 55+0 08 


1446.26+0.09 






72 


1,2,3,4 


0.41 


101.4 


591 27+0 08 


1438 58+0 08 


-19.0+5.8 


20.6+6.2 


73 


4 


0.41 


99.3 


883 47+0 08 


1332 41+0 09 






74 


4 


0.41 


95.0 


124.38+0.08 


477 45+0 08 






75 


1 


0.41 


127.7 


-108 53+0 08 


430 49+0 08 






76 


1,2 


0.41 


100.3 


12 37+0 08 


14 79+0 09 






77 


1 


0.41 


95.7 


62 90+0 08 

v./ — - . y vy _i_ v/ . v/ vy 


1271 66+0 08 






78 


1,2,3,4 


0.40 


101.8 


242 22+0 08 


1 240 44+0 1 


12 8+6 2 


18.9+7.3 


79 


1,2,3 


0.39 


99.6 


881 68+0 08 

U VJ X ilIU — i- V7 . VJ VJ 


1354 22+0 09 


1.5+11.6 


14.7+15.5 


80 


1,2,3,4 


0.37 


102.3 


623 10+0 08 

VJ.£._ > • J. y -L- \J . V/vj 


1 447 1 4+0 08 

X^^ / ■ X \7 . \7 O 


-7.7+6.1 


13.7+6.5 


81 


1 


0.37 


99.9 


871.98+0.08 


1307 43+0 10 






82 


4 


0.36 


93.0 


1 94 00+0 08 

X • \J \J _!_ \_/ . \J VJ 


476 56+0 09 






83 


1,2,4 


0.36 


100.7 


83.42+0.08 


274.56+0.10 


15.4+6.8 


-11.1+8.1 


84 


2,4 


0.35 


100.5 


822.73+0.08 


1426.15+0.10 






85 


1,2,3,4 


0.34 


96.7 


-623 32+0 08 


376 77+0 09 


-13.0+6.6 


0.0+7.5 


86 


4 


0.34 


99.5 


204 01+0 08 

£*\J^ 9 \J X _!_ \J . \J VJ 


1268 80+0 09 






87 


1 


0.34 


101.4 


-12 45+0 08 

X i_ . T ^ _!_ VJ . VJ VJ 


-3 26+0 08 






88 


4 


0.33 


100.0 


29 61+0 09 

7 . VJ X -l—\J.\Jy 


12 33+0 10 






89 


4 


0.33 


99.7 


203.21+0.09 


1285.90+0.11 






90 


4 


0.33 


98.6 


8.88+0.08 


4.77+0.09 






91 


1,2,3,4 


0.32 


104.9 


193.07+0.08 


782.37+0.09 


11.7+6.0 


-4.4+6.7 


92 


2,3,4 


0.31 


95.0 


-516.60+0.09 


256.92+0.09 


-12.9+9.6 


-12.7+10.4 


93 


3 


0.31 


94.6 


970.64+0.08 


1068.73+0.08 






94 


1,2,3 


0.31 


102.3 


859 58+0 08 


1412 65+0 09 

x r x .yj^s _!_ \j ,\j y 


-28.7+10.4 


51.8+12.4 


95 


4 


0.31 


99.1 


-85 10+0 08 

vj^j ^ x \J . \J vj 


96 72+0 09 






96 


1,2,3 


0.30 


123.6 


-91 29+0 08 

y X ■ ^ .7 _!_\J »\J VJ 


448.82+0.08 


-2.5+10.2 


-39.5+11.2 


97 


1,2,3,4 


0.30 


95.7 


5 1 fS 55+0 08 


265 90+0 08 

i-i\J*J • y \J -!—\J .\J yj 


-23.0+5.9 


3.6+7.0 


98 


3 


0.30 


99.3 


15.38+0.09 


8 48+0 09 






99 


2,3,4 


0.29 


95.1 


-616 35+0 08 

VJ X VJ * .7 _!_ \J . \J VJ 


433 17+0 09 

T^^^.X / -i-\J .\J y 


-12.4+9.3 


-9.3+10.3 


100 


1,2,3 


0.28 


95.8 


47 97+0 08 

~ / »J* r —1- VJ . VJ VJ 


194 52+0 08 

x y 1 . ^ ' — _!_ V7 . V7 ".7 


14.8+10.1 


-21.2+11.2 


101 


3 


0.28 


96.5 


970 57+0 08 

y 1 VJ 1 — !— UiVU 


1068 80+0 09 

i yj v7 \.> . v7 — i_ vy . vy y 






102 


4 


0.27 


90.6 


124.52+0.08 


477.19+0.09 






103" 


1,2,3 


0.27 


86.4 


106 46+0 08 


458 25+0 09 

1 ^ ' \ 1 . ^ — ./ _!_ V7 . V7 7 






104 


4 


0.27 


99.3 


-90 01+0 08 

yu.ui _!_ vy . vy vj 


142.78+0.10 






105 


4 


0.26 


98.6 


871 68+0 09 

yj i x > vj vj _!_ v/ . v/ ^ 


1282.10+0.10 






106 


1,3,4 


0.25 


101.5 


1 67+0 08 

1U.U / -J-VJ.VJVJ 


1208 55+0 09 

i — vy ^7 _!_ v / . vy y 


-32.1+6.3 


3.6+7.4 


107 


3 


0.24 


99.5 


203 91+0 09 

vj^j .y 1 J-v.u/ 


1268 48+0 10 

X V7 O . ^ Ky _!_ \7 • Iv 






108 


2 


0.24 


100.1 


-84.72+0.10 


97.24+0.12 






109 


2,3 


0.24 


95.7 


53 11+0 09 


1291 76+0 10 

x ^/ x ■ / vj _!_ v - x vy 






110 


1,2,3,4 


0.24 


96.4 


92 11+0 08 


467.69+0.09 


8.3+6.3 


-10.3+7.2 


111 


1,2 


0.23 


102.4 


-94 09+0 08 

_V^.VJ7_!_VJ.vJU 


164 25+0 08 

X V7T7Z_».' _!_V7 .\7 (J 






112" 


1,2,3 


0.22 


87.4 


41 52+0 08 

I X * .7 — _!_ V/ * V/ VJ 


525 37+0 09 

*7^%7-%7 ' — !_ V / . V / 7 






113 


2 


0.22 


100.5 


596 04+0 10 

JyUt vj . x vj 


1485.14+0.14 






114 


4 


0.21 


115.7 


47 ig+o 08 

I / . X VJ _!_ UiVU 


1065 65+0 09 

X \_/V_7.7 . V ' — !_ V7 . V./ 7 






115 


2 


0.21 


99.3 


11.83+0.11 


106 83+0 10 

X Vy Vy . VJ ^ _!_ vy . X V7 






116 


1 


0.20 


101.7 


gg 99+0 08 

y \j * y y _!_ vj .vj vj 


1384 97+0 09 

x ^y u ~ . y / _!_ yj *yj y 






117 


2,3 


0.20 


95.1 


1 35 26+0 09 

X .7 .7 . ^ VJ _!_ V./ . VJ V 


558 19+0 09 

.7.jtj. iyj-u.uy 






118 


4 


0.20 


95.9 


123.38+0.08 


475.38+0.09 






119 


3,4 


0.19 


95.7 


121.26+0.08 


475.30+0.09 






120 


2,3,4 


0.19 


95.7 


51.23+0.08 


1291 83+0 09 

x w y x ■ vy ^y _!_ yj *yj y 


-23.9+10.5 


52.2+12.1 


121 


4 


0.18 


91.7 


1 10 77+0 08 

X X VJ . / / _!_ \J . \J VJ 


460 47+0 09 

^\J\J / _!_ vy . V7 v 






122 


1,2,3,4 


0.18 


103.6 


-54 77+0 08 


1049.45+0.09 


-16.3+6.5 


14.6+8.1 


123 


1,2,3,4 


0.17 


94.3 


105 02+0 08 

X VJ^/ . VJ _!_ V/ . \J O 


1351 io+0 10 

X X . X V/_!_V/. X \J 


-18.2+6.5 


15.0+7.7 


124 


1,2 


0.17 


116.3 


ie 99+0 08 

L \j • y y -i-\J .\J \J 


1 1 62 29+0 09 

X X VJ . 2— y _!_ \J ,\J y 






125 


2 


0.16 


101.2 


-17.73+0.14 


-31.34+0.15 






126 


1 


0.16 


92.3 


120.38+0.08 


469.82+0.09 






127 


1 


0.16 


101.6 


-53.24+0.08 


1046.79+0.09 






128 


1,3 


0.16 


101.7 


194.02+0.09 


1216.44+0.10 






129 


1 


0.15 


91.0 


-227.30+0.08 


646.92+0.10 






130 


3,4 


0.14 


91.7 


-167.73+0.09 


380.68+0.10 






131 


4 


0.14 


84.2 


-507.22+0.08 


520.81+0.10 







20 



L. Moscadelli et al.: A Double- Jet System in the G3 1.4 1+0.31 Hot Molecular Core. 

Table 4. continued. 



Feature 


Epochs of 


I k 


Vlsr 


A x 


A v 


v x 


V y 


Number 


Detection 


(Jy beam ) 


(km s~') 


(mas) 


(mas) 


(km s~') 


(km s" 1 ) 


132 


4 


0.14 


111.2 


20.79±0.08 


1155.17+0.09 






133 


2 


0.14 


106.7 


-105 08+0 09 


456 25+0 10 






134 


1,2,3 


0.14 


105.0 


-407 55+0 08 


89 72+0 10 


-89.4+11.3 


-60 6+13 6 

VJVJ • \7— l_ X *7 ■ V7 


135 


1 


0.14 


88.7 


123 51+0 08 


4S7 91+0 10 






136 


1 


0.14 


94.0 


_7 07+0 08 

/ . VJ I _!_V/.V/U 


1 203 64+0 09 






137 


2 


0.14 


101.9 


-17.06+0.10 


_30 70+0 10 

*7 \J ■ / \7 — !— V7 » ± V 






138 


3 


0.13 


115.0 


42.78+0.08 


1076.47+0.09 






139 


1,2 


0.13 


104.5 


-52 94+0 08 


1062.76+0.10 






140 


3 


0.13 


87.3 


108.27+0.08 


506 54+0 09 






141 


2 


0.13 


96.5 


131.24+0.10 


488 19+0 11 

"UUt x y _!_ vj ■ x x 






142 


1,2 


0.13 


91.7 


-92 49+0 08 


677 89+0 09 

\j i i *\j y _!_ V7 • \7 v 






143 


1,2,3,4 


0.13 


87.0 


-62 91+0 08 


279 80+0 10 

j-j t y . vj vj — !— V7 tiv 


22.6+8.0 


-62.4+9.7 


144 


1,2,3,4 


0.12 


104.7 


-4S 49+0 08 

~-J J -i-\J .\J \j 


1082 40+0 09 


-27.6+7.4 


-2.5+9.6 


145 


2 


0.12 


118.5 


35 40+0 09 


1184 50+0 10 






146 


2,3 


0.11 


124.0 


17 42+0 09 

X / . ^ _!_ \J » \J y 


1 158 01+0 11 

X 1 JUtVl _!_ \J . X X 






147 


4 


0.11 


90.1 


103 48+0 09 


458.33+0.11 






148 


3,4 


0.11 


103.4 


33 25+0 10 

mj *J * £*^J _I_\7 » X V/ 


1245.23+0.11 






149 


2 


0.11 


93.7 


-12 45+0 09 

l j— u r ^ _!_ vj .vj y 


-3.75+0.11 






150 


2 


0.11 


88.6 


106 Q^+O 09 


458 31+0 10 






151 


2 


0.11 


117.4 


50 55+0 10 

.7 \7 . ,7 .7 _!_ V.lu 


H90 01+0 10 






152 


1 


0.11 


96.3 


_0 54+0 09 

\j . ^> ^ j-vj .vj y 


1 199 17+0 10 

xx.y.y.x /_!_vy.xv/ 






153 


1 


0.10 


106.4 


_95 93+0 08 

y ^' . y *j _!_ v.' . v.' 


163 34+0 10 

1 V J - .7 r_!_ VJ* 1 V 






154 


2 


0.10 


102.0 


-53 00+0 10 


1046 06+0 13 

X \7 I V7 • \J\J -1— V.J . X .7 






155 


1,2,3,4 


0.10 


104.3 


30 90+0 08 

JVi 7 V./ _!_ V./ . V / V.I 


1237.90+0.09 


-14.9+7.6 


16.7+9.2 


156 


1,2,3 


0.10 


92.1 


-93.50+0.09 


174.90+0.12 


-23.3+12.0 


-7.4+15.1 


157 


1,2 


0.09 


93.7 


-109.34+0.08 


170.88+0.10 






158 


2 


0.09 


102.6 


119.90+0.10 


475.82+0.11 






159 


4 


0.09 


115.4 


22.03+0.10 


1155.46+0.12 






160 


1 


0.09 


103.5 


597 29+0 09 


987 29+0 1 1 

y vj i . y _l_ vy# i i 






161 


1,2 


0.09 


104.8 


_56 70+0 09 

^i vj . / vj-<-\j *\j y 


1041 70+0 12 

X VJT X . / \J ■ X ^ 






162 


1,2 


0.09 


92.7 


167 30+0 08 

X \J 1 * \J _!_ \7 . \7 VJ 


1841.14+0.10 






163 


2,4 


0.09 


87.9 


106 06+0 10 

x vjwj. \j\j-L- \j. x \j 


457 91+0 11 

^»J / . y X _!_V7 . X X 






164 


2 


0.08 


92.3 


104.04+0.11 


458.34+0.12 






165 


1 


0.08 


112.0 


-106 06+0 08 

1UU . \7 V7 _!_ 17 » 17 VJ 


409 95+0 1 1 

~vj y . y ~j _!_ v/ • x x 






166 


1,2 


0.08 


91.6 


151.34+0.08 


1831.39+0.11 






167 


3 


0.07 


102.9 


25.19+0.12 


1241.43+0.12 






168 


1 


0.07 


92.9 


-90.08+0.09 


678.15+0.11 






169 


1 


0.06 


113.2 


-0.62+0.09 


1063.01+0.11 






170 


1 


0.06 


93.4 


139.41+0.09 


518.74+0.11 






171 


1 


0.05 


111.8 


-106.28+0.11 


390.79+0.13 






172 


1 


0.05 


88.7 


103.29+0.11 


458.90+0.12 






173 


2 


0.05 


113.8 


3.78+0.12 


1128.38+0.15 







Notes. Column 1 gives the feature label number; column 2 lists the observing epochs at which the feature was detected; columns 3 and 4 provide 
the intensity of the strongest spot and the intensity-weighted LSR velocity, respectively, averaged over the observing epochs; columns 5 and 6 give 
the position offsets (with the associated errors) along the R.A. and Dec. axes, relative to the feature #1, measured at the first epoch of detection; 
columns 7 and 8 give the components of the relative proper motion (with the associated errors) along the R.A. and Dec. axes, measured with 
respect to the reference feature #0 (the "center of motion"). 
{a) Feature with not reliable proper motion, even if observed at three epochs. 
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Table 5. 6.7 GHz CH 3 OH Maser Parameters 



Feature 


Epochs of 


Ipeak 


^LSR 


A x 


Ay 




Vy 


Number 


Detection 


(Jy beam -1 ) 


(km s ) 


(mas) 


(mas) 


(km s ) 


n -1 \ 

(km s ) 





1,2,3,4 




97.5 


192.50±0.06 


-399.23+0.07 


0.0+0.0 


0.0+0.0 


1 


1,2,3,4 


9.57 


92.3 


0.00±0.00 


0.00+0.00 


2.8+1.1 


3.4+1.2 


2" 


1,2,3,4 


2.64 


104.2 


347.78±0.08 


-47.91+0.09 






3 


1,2,3,4 


2.48 


103.4 


346.9 1±0.08 


-38.57+0.08 


-5.9+1.0 


15.6+1.2 


4 


1,2,3,4 


1.66 


95.9 


-241.15+0.09 


-85.27+0.09 


18.2+1.3 


-3.0+1.3 


5 


1,2,3,4 


1.63 


103.6 


356.28±0.08 


-71.34+0.08 


-18.9+1.2 


18.7+1.4 


6 


1,2,3,4 


1.14 


93.8 


25.76±0.08 


-886.78+0.09 


0.5+1.1 


-8.1 + 1.3 


1" 


1,2,3,4 


0.96 


95.6 


-254.83±0.09 


-75.53+0.09 






8 


1,2,3,4 


0.95 


102.9 


357.14±0.08 


-65.98+0.09 


-8.9+1.1 


16.9+1.4 


9 


1,2,3,4 


0.86 


96.0 


596.01±0.09 


-704.83+0.10 


11.1+1.3 


-11.5+1.6 


10 


1,2,3,4 


0.84 


105.4 


340.86±0.08 


-4.30+0.09 


-2.2+1.4 


10.7+1.7 


11 


1,2,3,4 


0.80 


95.6 


591.89±0.09 


-674.21+0.10 


11.3+1.3 


-10.1 + 1.5 


12 


1,2,3,4 


0.77 


102.7 


353.48+0.08 


-51.70+0.09 


-7.3+1.2 


14.8+1.5 


13 


1,2,3,4 


0.75 


96.3 


543.18+0.09 


-1092.66+0.10 


4.9+1.3 


-7.6+1.6 


14 


1 


0.75 


103.5 


349.46+0.10 


-54.71+0.11 






15 


1,2,3,4 


0.74 


95.2 


-260.43+0.10 


-72.04+0.10 


14.5+1.3 


-4.6+1.4 


16 


1,2,3,4 


0.71 


96.3 


-134.94+0.11 


-120.17+0.11 


-4.8+1.8 


6.2+1.9 


17 


1,2,3,4 


0.61 


98.8 


624.01+0.09 


-731.91+0.10 


4.9+1.4 


-4.7+1.7 


18" 


1,2,3,4 


0.52 


99.1 


177.27+0.10 


136.98+0.13 






19 


1,2,3,4 


0.42 


99.2 


171.80+0.12 


147.41+0.14 


1.5+1.7 


6.0+2.2 


20 


1,2,3,4 


0.35 


93.3 


55.71+0.11 


-691.54+0.12 


0.7+1.6 


-9.9+2.0 


21 


1,2,3,4 


0.32 


96.8 


-148.45+0.13 


-114.68+0.13 


-0.8+1.9 


4.9+2.2 


22" 


1,2,3 


0.31 


98.7 


616.63+0.10 


-729.16+0.11 






23 


1,2,3,4 


0.29 


97.6 


604.58+0.12 


-711.77+0.13 


-1.9+1.8 


-3.3+2.2 


24 


1,2,3,4 


0.24 


106.6 


334.10+0.14 


-20.74+0.18 


-13.8+1.8 


17.8+2.5 


25 


1,2,3,4 


0.23 


97.2 


75.48+0.18 


181.85+0.29 


1.5+2.4 


4.4+4.2 


26 


1,2,3,4 


0.23 


105.1 


339.26+0.11 


-68.06+0.13 


-6.6+2.1 


17.8+2.7 


27 


1,2,3,4 


0.20 


92.6 


91.99+0.16 


-683.14+0.18 


7.6+2.5 


-15.4+3.1 


28 


1,2 


0.19 


105.5 


338.14+0.15 


-70.04+0.19 






29 


1,2 


0.18 


96.1 


-105.40+0.19 


-151.75+0.22 






30 


1,2,3,4 


0.17 


98.2 


118.15+0.24 


116.94+0.38 


3.3+3.1 


3.8+4.7 


31 


4 


0.17 


98.5 


614.03+0.18 


-728.71+0.23 






32 


1,2 


0.17 


96.2 


-263.36+0.24 


-78.72+0.21 






33 


1,2,3 


0.16 


98.4 


612.57+0.16 


-727.02+0.17 


-1.9+5.0 


-8.6+5.8 


34 


1,2,3,4 


0.16 


94.9 


-1.41+0.18 


-989.38+0.17 


-1.9+3.0 


-5.9+3.3 


35 


1,2,4 


0.16 


95.4 


577.99+0.27 


-669.69+0.32 


21.0+4.2 


-10.1+4.8 


36" 


1,2,3,4 


0.16 


94.4 


27.33+0.20 


-990.87+0.26 






37 


1,2,3,4 


0.15 


97.3 


-133.29+0.24 


-128.33+0.27 


-6.7+3.4 


9.0+4.0 


38 


1,2,3,4 


0.15 


92.6 


57.45+0.21 


-684.51+0.22 


1.8+3.2 


-14.1+3.7 


39" 


1,2,3 


0.15 


94.6 


22.81+0.16 


-983.89+0.20 






40 


1,2,3,4 


0.15 


99.8 


635.73+0.19 


-746.47+0.25 


1.3+3.1 


-4.1+4.0 


41 


1,2,3,4 


0.14 


100.0 


714.67+0.17 


-1004.55+0.19 


5.4+3.9 


-8.9+4.2 


42" 


1,2,3,4 


0.14 


93.0 


36.56+0.30 


-694.89+0.28 






43 


1,2,3 


0.13 


95.2 


14.57+0.25 


-978.73+0.23 


-0.4+8.2 


-22.1+8.6 


44 


1,2 


0.13 


95.9 


-95.07+0.28 


-138.93+0.37 






45 


1 


0.13 


95.8 


587.12+0.28 


-667.66+0.54 






46 


1,2,3,4 


0.11 


94.7 


-341.89+0.16 


-474.74+0.20 


-5.4+4.0 


0.0+5.1 


47 


1,2,3,4 


0.10 


94.3 


-113.82+0.22 


-489.24+0.28 


-0.3+3.7 


-9.6+5.0 


48 


1,2,3,4 


0.10 


94.5 


14.55+0.17 


-989.79+0.20 


-11.1+4.6 


-8.0+5.4 


49 


1,2,3,4 


0.10 


92.0 


57.76+0.25 


-825.33+0.36 


5.9+4.5 


-14.3+6.5 


50 


2 


0.09 


93.8 


18.28+0.25 


41.24+0.41 






51 


3 


0.09 


98.4 


104.20+0.29 


137.52+0.43 






52 


2 


0.09 


94.5 


14.85+0.29 


-1051.79+0.35 






53" 


1,2,3,4 


0.08 


95.0 


-270.45+0.40 


-71.69+0.31 






54 


1,2,3,4 


0.08 


93.0 


47.88+0.18 


-690.79+0.22 


1.2+5.1 


-10.7+6.0 


55" 


1,2,3,4 


0.08 


94.3 


-9.62+0.39 


100.93+0.62 






56 


1,2,3,4 


0.08 


94.9 


572.04+0.31 


-658.57+0.43 


4.3+5.2 


-8.7+6.6 


57" 


1,2,3,4 


0.08 


106.7 


320.07+0.19 


-86.91+0.28 






58 


1 


0.08 


94.9 


8.04+0.39 


-988.51+0.47 






59 


1,2 


0.08 


93.9 


134.44+0.23 


-730.51+0.35 






60 


1,2,3,4 


0.08 


97.9 


86.04+0.38 


152.20+0.57 


2.2+5.1 


11.6+8.0 


61 


1,2,3 


0.07 


101.0 


727.27+0.31 


-917.53+0.38 


13.4+9.0 


-19.4+11.1 


62 


1,2,3 


0.07 


97.0 


53.64+0.36 


177.40+0.59 


1.2+8.2 


5.0+13.2 


63 


1,2 


0.07 


104.5 


301.71+0.25 


-80.64+0.42 






64 


1 


0.07 


96.5 


-48.64+0.35 


-844.58+0.48 






65 


1 


0.06 


94.5 


-122.40+0.39 


-492.51+0.43 
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Table 5. continued. 



Feature 


Epochs of 


I k 


Vlsr 


A x 


A v 




Number 


Detection 


(Jy beam ^ ) 


(km s~') 


(mas) 


(mas) 


(km s (km s 


66 


1 


0.06 


93.2 


604 31+0 37 


-664 64+0 48 




67 


2 


0.06 


98.6 


103 39+0 62 


142.21+1.28 




68 


4 


0.06 


98.5 


103 72+0 76 


139 63+r 98 




69 


1 


0.06 


95.0 


98 92+0 40 


-822.30+0.52 




70" 


1,3,4 


0.06 


97.0 


-239 70+0 49 


-82 07+0 68 




71 


1 


0.05 


98.6 


102.52+0.47 


144.12+0.88 




72 


1,2 


0.05 


93.6 


-17.51+0.39 


38 20+0 64 




73 


1 


0.05 


93.8 


115 21+0 35 


-1317 51+0 56 




74 


1 


0.05 


106.6 


340.24+0.29 


-5 02+0 45 




75 


1,2 


0.05 


93.1 


98 92+0 69 


-676 75+0 64 




76 


1,2 


0.05 


93.6 


-8 51+0 39 


13 39+0 72 




11" 


1,2,4 


0.05 


98.0 


527 33+0 40 


-726 68+0 66 




78 


' 1 


0.05 


94.2 


-6.38+0.46 


204.54+0.63 




79 


2 


0.05 


91.5 


33.46+0.38 


56.86+0.66 




80 


1 


0.05 


92.7 


83.34+0.59 


-690.23+0.80 




81 


3 


0.05 


105.8 


351.51+0.53 


-160.45+0.67 




82 


1 


0.04 


94.3 


-9.02+0.48 


213.87+0.73 




83 


3 


0.04 


94.7 


-46.01+0.71 


146.18+0.84 




84 


1 


0.04 


90.5 


-137.65+0.46 


-532.52+0.77 




85 


1 


0.04 


97.0 


-239.12+0.66 


-91.77+0.63 





Notes. Column 1 gives the feature label number; column 2 lists the observing epochs at which the feature was detected; columns 3 and 4 provide 
the intensity of the strongest spot and the intensity-weighted LSR velocity, respectively, averaged over the observing epochs; columns 5 and 6 give 
the position offsets (with the associated errors) along the R.A. and Dec. axes, relative to the feature #1, measured at the first epoch of detection; 
columns 7 and 8 give the components of the relative proper motion (with the associated errors) along the R.A. and Dec. axes, measured with 
respect to the reference feature #0 (the "center of motion"). 

{a) Feature with not reliable proper motion, even if observed at three or four epochs. 
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